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SCOUT TRAJECTORY ERROR PROPAGATION 
COMPUTER PROGRAM 


by T, R. ^fyle^ 
Vought Corporation 


■ ^ SUMMARY 

This report describes a FORTRAN coded computer program which calculates 
trajectory error covariance matrices and statistics from a data base of flight 
experience* The data base consists of trajectory errors resulting from past 
flights. A covariance matrix is calculated and may be propagated in time and 
added to a boost covariance matrix of a spin-stabilized stage* A sensitivity 
matrix is developed from the eigenvalues and eigenvectors of the final 
propagated covariance matrix. This sensitivity matrix is provided for use by 
another computer program to randomly sample the matrix using a Monte Carlo 
technique to yield sample errors which will produce the same covariance matrix. 

The theory and methods presented in this report for calculating error 
statistics and propagating an error covariance matrix are of general interest 
since they have applications other than those contained herein* 

Included in this report are program theory, user instructions, output 
descriptions, subroutine descriptions and detailed FORTRAN coding information* 


- 1 - 



1.0 lOTRODOCTION 


Since 1969, flight experience has been used as the basis for predicting 
Scout orbital accuracy. _ The data base used for calculating the accuracy 
consists of errors in the trajectory parameters (altitude, velocity, etc.) at 
stage burnout as observed on Scout flights. Approximately 50 sets of errors 
are used in a Monte Carlo analysis to generate error statistics in the 
trajectory parameters. A covariance matrix is formed which may be propagated 
in time. The mechanization of this process resulted in computer program Scout 
Trajectory Error Propagation (acronym STEP) and is described herein. 

Computer program STEP may be used in conjunction with the Statistical 
Orbital Analysis Routine (Reference 1) to generate accuracy in the orbit 
parameters (apogee, perigee, inclination, etc.) based upon flight experience. 



2.0 DEFINITIONS 


2.1 NoCatioa 

Symbols used in this report are listed below with their definitioa and 
units. 

English Alphabet 
Ce 

Cn 

CR 
h 
n 

^e 
R 

V 


Flight path angle, deg 
Deviation from nominal 
Random number, unitless 
Velocity azimuth, deg 
Latitude, deg 
Longitude, deg 
Statistical mean 
Correlation coefficient 
Summation 
Standard deviation 


Adjusted value 

Partial derivative of x to y 
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Greek Alphabet 

y 

A 

e 

r 

X 

d 

w 

p 

2. 

a 

Others 

/ 

3X 

*57 


East component of position error, ft 
North component of position error, ft 
Crossrange, n.mi. 

Altitude, n.mi. 

Number of samples 
Earth radius, ft 
Range, n.mi. 

Velocity, fps 



2.2 Flowchart Conventions 


Flowchart conventions used in this report are as follows 


Process 


Input/Output 




Svibroutine 
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3.0 PROGRAM DESCRIPTION 


This section describes the utilization of the Scout data base, program 
theory, user instructions and output definitions. 

3.1 General 

The putiJose of computer program STEP is to calculate trajectory error 
statistics from a data base of errors resulting from actual flights. This 
data base consists of stage burnout errors in the trajectory parameters - 
altitude, velocity, flight path angle, flight azimuth, latitude and longitude. 

The Scout data base at orbit insertion is shown in Table 3.1. These 
errors represent flights at various launch azimuths and insertion altitudes. 
Since errors in latitude and longitude are dependent upon launch azimuth and 
since errors in altitude, velocity and flight path angle may be dependent upon 
insertion altitude, adjustments are made to the flight errors in order to 
provide a consistent data base from which to generate the error statistics. 
These adjustments are discussed in the following paragraphs. 

The Scout data base includes flights with launch azimuths ranging from 
easterly to southerly to slightly west of south. Latitude and longitude 
errors, which may produce inclination errors, result from crossrange errors 
on easterly flights and range errors on southerly flights. Therefore, 
latitude and longitude errors of the flight history data base are not a 
consistent set of sample errors for a given launch azimuth. Since it 
Is necessary to have a consistent set of sample errors, Che data base latitude 
and longitude errors are adjusted to the launch azimuth of interest. From 
the nominal flight azimuth, latitude, longitude and their errors on each 
flight of the data base, range And crossrange errors are calculated. Since 
range and crossrange errors are independent of the flight azimuth, these 
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errors caa be converted to new latitude and longitude errors for the flight 
azimuth of interest by the following transformation. 

Cj^= fgAX 
C£= fgAycosX 
AR “ f + Cg sin f 

ACR= sm f - Cg cos f 

The' above relationships are evaluated for each flight sample of the data 
base, resulting in range and crossrange errors for each flight. Note that 
latitude and longitude errors are needed to calculate range and crossrange 
errors. The range and crossrange errors are then converted to latitude and 
longitude errors applicable to the flight azimuth of interest as shown below. 

* A R cos f '+ A CR sm 
Cg = A R sm f' — A CR cos f ' 

AX'=CN7re 
A n' = C^'lijcos X 


where the primed values pertain to the conditions at which the error 
statistics are desired. 

The above process yields latitude and longitude errors which are used in 
the flight data base for calculating trajectory error statistics. 

The second adjustment made to the flight data base is to "normalize" the 
errors in altitude, velocity and flight path angle to the insertion altitude 
of interest. The purpose of this adjustment is to obtain the flight errors 
which would have resulted if all the flights had had the same insertion 
altitude. This adjustment is accomplished for each flight sample as follows; 


AV'= — 
ov 
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AV 



where AV'« adjusted velocity error at the insertion altitude of interest 
AV» flight sample velocity error 

09' * standard deviation of velocity at the insertion altitude of 
interest 

* standard deviation of velocity at the flight sample insertion 
altitude 

Similar expressions are used for altitude and flight path angle errors. 
Azimuth, latitude and longitude errors are not adjusted for insertion altitude 
because they are independent of insertion altitude. 

Ihe "normalized" deviations in altitude, velocity and flight path angle -* 
the adjusted deviations in latitude and longitude - and the azimuth deviations 
from the flight samples - are combined to be the flight data base from which 
error statistics are calculated. This flight data base is applicable only to 
a trajectory with the insertion altitude and launch azimuth to which the 
flight errors were adjusted. 

The "normalizing" process should be used only when altitude, velocity and 
flight path angle errors input in the flight data base, are a function of 
insertion altitude. This relationship is true for Scout because there is a 
long coast time (300-400 seconds) prior to the last stage boost, which allows 
the errors to grow. Thus, the error magnitudes at last stage ignition are a 
function of time and altitude. If the input flight data base corresponds to 
the stage burnout prior to the long coast time, the errors should not be 
normalized. 

The normalizing equation involves a ratio of standard deviations of the 
trajectory parameter. These deviations can be obtained by inputting a flight 
data base at stage burnout prior to the long coast time and propagating the 
covariance matrix to various insertion altitudes. 
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3.2 Program Theory 


The primary function of STEP is Co calculate and/or propagate a covariance 
matrix of trajectory state parameters - altitude, velocity, flight path angle, 
azimuth, latitude and longitude. The initial covariance matrix is either 
input directly or is calculated from input samples of actual flight errors 
(flight experience). Covariance propagation is optional and is controlled by 
inputs to the program. If selected, a boost covariance matrix of a 
spin-stabilized stage is calculated and added to the propagated covariance 
matrix. The resulting matrix - or the covariance matrix at the input epoch if 
it is not propagated - provides statistics in the trajectory parameter 
errors. Also, sensitivity coefficients calculated from the covariance matrix 
can be used by another computer program to provide a sampling of trajectory 
errors which can be converted to orbital parameter error statistics. Such a 
program is the Statistical Orbital Analysis Routine, Reference (I). 

The calculation of a covariance matrix from flight results and the 
calculation of the sensitivity coefficients are discussed below. 

Each of the six parameters of the flight data base, as described in 
Section 3.1, are combined in the manner shown below to yield mean values, 
standard deviations and correlation coefficients. 


S" Axi 
a. 



naxey 


mean value 


standard deviation 


correlation coefficient 


where x and y are any two parameters 
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The error covariaace matrix is obtained from the error statistics 


calculated as shown on the preceding page and has the following form: 


0}j2 Synmetric 

Phv°h*^v 

PhK^: Pv^v^ 

PhX°h®X PvX^yOX 
f^hn%9n PsyP'fin 


The error covariance matrix is a real, symmetric matrix and can be 


diagonalized to obtain the eigenvectors and eigenvalues. A matrix using the 


eigenvalues is formed as follows: 



where " eigenvalues 

The sensitivity matrix is formed from the A matrix and the eigenvector 
matrix as follows: 


[S] = [AHET]T 

where is the transpose of the eigenvector matrix. 

The S matrix represents a sensitivity matrix of the six trajectory 
parameters to six independent and uncorrelated error sources and has the 

following form: 

9h dV dj dj d\ du 

dE| 9E| dEj 9E^ 9Ej 9E| 

9h 9V dy dj 9\ du 

9E2 9E2 9E2 9E2 9E2 9E2 

9h 9V dy dj d\ du 

9E3 9E3 9E3 9E3 9E3 9E2 

[S] = 

9E4 9E4 9E4 9E4 9E4 9E4 

9h 9V dy gf 9X 9^ 

9E5 9E5 9E5 9E5 9Ej 9E5 

9h 9V 97 ^ ^ ^ 

9E^ ^ 9Efi 9Efi 9Eg 





PXji^X'V 
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where » Independent, uncorrelated error sources 
The sensitivity matrix is used to calculate random errors in the 
trajectory parameters as follows: 


Ah 

AV 

Ay 

Ar 

AX 

Ay 



where are random numbers from a normal distribution with a 
mean value of zero and a standard deviation of one 


A covariance matrix is propagated by using a Monte Carlo analysis to: (1) 

sample the trajectory errors from the sensitivity matrix, (2) add the errors 
to the nominal values, (3) propagate the state parameters along a conic by a 
time Increment, and (4) subtract the new state parameters from the nominal 
values at the new epoch* The resulting trajectory errors are combined as 
described on the preceding pages to yield the mean values, standard deviations 
and correlation coefficients and, thus, the covariance matrix* 

The spin-stabilized stage boost covariance matrix is also formed using 
a Monte Carlo analysis* The boost error sources are sampled as shown below 
and combined as shown previously in this section* 


Ah 


9h 

3h 

3h 

3h“ 

€1 


3 Ei 

3 E 2 

3E3 

3E4 


AV 


3V 


3V 

3V 




3 Ei 

3E2 

3E3 

3E4 


At 


37 

37 

37 

37 



3^^ 

3E2 


3E4 


Ar 


3 r 

3 r 

3 r 

3 r 



3 Ej 

^2 

3F3 

3E4 


AX 


3X 

3X 

3X 

3X 



3 Ej 

3E2 

3E3 

3E4 


Ay 


3y 

3y 

3y 

3y 



3E^ 

W2 

3E3 

3E4 

^4 
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when E| 

i2L 

aEj 


X 


"i 


= boost error sources 

= partial derivative of respective 
trajectory parameter to a one sigma 
magmtude of the error source 

= a random number defined as follows 
For motor performance error source 
^l=Xnj 

For pitch tipoff error source 

«2-IXn2l^o*(Xu*360“) 

e3»lXn2lan(Xu*360*) 

«4 = Xn3 

= a normally distnbuted random number 
with a mean value of 0 and a standard 
deviation of 1. 


Xy = a uniformly distributed random 
number from 0 to 1. 


3 » 3 User Instructions 

STEP utilizes both fixed field and a modified FORTRAN NAMELIST for data 
input. The flight experience data base is input via fixed field format since 
the data base normally does not change* The remaining data, which normally 
changes with each program execution, is input via NAMELIST in order to utilize 
a flexible input format. 

3.3.1 Fixed Field Format - The flight data base is input via fixed field 
format. For each flight the following parruaeters are input at an epoch, 
o flight identification 

o option specifying usage of latitude, longitude euid azimuth errors* 
o nominal altitude 
o nominal latitude 
o nominal flight azimuth 

o observed errors in altitude, velocity, flight path angle, flight 
azimuth, latitude, longitude. 
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Fixed field data must either be right justified in the field or contain a 
decimal point. 

The flight experience data base is input in the order shown below 
for each flight. Data may be input for up to 100 flights. Preceding the 
flight data, one title card is used for identification purposes only. 
Following the flight data, an END card beginning in column 1 is used to 
specify the end of the flight data. 


Column Range 
1-5 
6-9 


10-16 

17-23 

24-30 

31-37 

38-44 

45-52 

53-60 

61-70 

71-80 


Fixed Field Input Definition 

Definition 
Flight Number 
Flag 

* 1 if flight is to be used for latitude and 
longitude statistics only. 

^ 2 if flight is to be used for all statistics. 

» 3 if flight is to be used for azimuth statistics 
only. 

4 is flight is not to be used for statistics. 
Nominal altitude, n.mi. 

Altitude error, n.mi. 

Velocity error, fps 
Path angle error, deg 
Azimuth error, deg 
Latitude error, deg 
Longitude error, deg 
Nominal latitude, deg 
Nominal azimuth, deg 
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3.3.2 NAMELIST Format - A modified FORTRAN NAMELIST is used for inputting 


data to STEP. NAMELIST is used because of its readability and simplicity of 
inputting data. The following rules apply to NAMELIST input to STEP. 

1. First card of a data group or case is (INPUTD beginning in column 2. 
Blanks are not allowed. 

2. Last card of a data group or case is ^END beginning in column 2. 
Blanks are not allowed. 

3. Blanks may not be used within names but may be used elsewhere. 

4. Variable names are followed by an equal sign, followed by a value, 
followed by a comma, e.g., NSAMP^IOOO, 

5. Only columns 2*-72, inclusive, are used. 

6. Titling information may be input by the appropriate title names, 
e.g., TITLEl- LOW ALTITUDE TRAJECTORY ERROR STATISTICS 

TITLEl must begin in column 2. 

7. Any number of names and values may be on a single card or line. 

8. Complete data arrays are input in the following form: 
name * value, value, value, ..., 

Data values may be continued on the next line, but the last 
character on every line must be a comma, excluding title cards. 

9. Repeated data values may be input by using a repetition factor and 
an asterisk, e.g., DATAG • 0.4, 4*0.45, 0.5, 0.65, 

10. One or more specific elements of an array may be input, e.g., 

EMAG(2)- 1.2, 1.6, 

Subsequent data cases are allowed by providing additional sets of NAMELIST 
data. All input data is retained for subsequent cases but can be changed by 
inputting new values. 

A sample data case is included as Appendix A to exemplify data case setup. 
Definitions of specific NAMELIST inputs to STEP are shown below. Default 
values are shown when they are set by the program prior to reading input 
data. Data units are feet, degrees and seconds unless otherwise noted. 
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NAMELIST Input Definitions 


COV Covariance matrix of altitude, velocity, flight 

path angle, azimuth, latitude and longitude. 

Input when lERROR =» 4 

Standard deviation in flight path angle, used to 
normalize flight errors. Array of 7 values. 
Independent variable is HINJ. 

(0.428, 0.428, 0.428, 0.430, 0.434, 0.442, 0.453 
built-in) 

Standard deviation in altitude, used to normalize 
flight errors. Array of 7 values. Independent 
variable is HINJ. 

(3.35, 4.85, 5.80, 6.45, 6.90, 7.20, 7.30 
built-in) 

Standard deviation in velocity, used to normalize 
flight errors. Array of 7 values. 

Inde pendent 18ariable is HINJ. 

(75.6, 84.6, 88.8, 91.3, 92.5, 93.0, 93.0 
built-in) 

Values for error sources used to calculate 
covariance matrix of a spin stabilized stage. 
Array of 4 values. Input when propagation is 
used. 

(1) * ratio of standard deviation desired to the 
standard deviation used to determine SENl 
for the motor performance error source. 

(2) > same as above except for determining SEN2 
for the pitch tipoff error source. 

(3) * same as above except for determining SEN3 
for the yaw tipoff error source. 

(4) same as above except for determining SEN4 
for the timer error source. 

HINJ Altitude used to normalize flight errors. Array 

of 7 values. Dependent variables are DATAG, 
DATAH, - DATAV. 

(100., 200., 300., 400., 500., 600., 700. 
built-in) 

lERROR Option for inputting data errors 

■ 1 Input flight results of altitude, velocity 
and path angle errors, 100 samples or less. 
Errors are input in fixed field format prior 
to NAMELIST data. 

“ 2 Input flight results of altitude, velocity 
path angle, azimuth, latitude and longitude, 
100 samples or less. Errors are input in 
fixed field format prior to NAMELIST data. 

(2 built-in) 


DATAG 


DATAH 


DATAV 


EMAG 
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NAMELIST Input Definitions (Continued) 

* 3 Input standard deviations (SIG) and 

correlation coefficients (RHO) of altitude, 
velocity, path angle azimuth, latitude and 
longitude. 

■ 4 Input covariance matrix (COV) of altlttxde, 
velocity, path angle, azimuth, latitude and 
longitude. 

NERBOR 

Number of error sources of the spin stabilized 
stage. (3 built-in) 

NORM 

Non-zero value normalizes altitude, velocity and 
path angle errors of the flight data base to the 
altitude of SI. (1 built-in) 

NSAMF 

Number of samples used In Monte Carlo analyses. 
(5000 built-in) 

RHO 

Correlation coefficients In order of altitude 
velocity, path angle, azimuth, latitude and 
longitude. Array of 15 values. Input when 
lERROR - 3. 

SENl 

Sensitivity of spin stabilized stage burnout 
state parameters to one sigma motor performance 
error source. Units are state parameter 
unlts/slgma. Array of 6 values. Order of state 
parameters are altitude, velocity, path angle, 
azimuth, latitude and longitude. 

m 

SEN2 

Same as SENl except error source Is pitch 
tlpoff. Units are state parameter unlts/deg. 

SEN3 

Same as SENl except error source Is yaw tlpoff. 
Units are state parameter unlts/deg. 

SEN4 

Same as SENl except error source Is stage 
Ignition time. Units are state parameter 
unlts/sec. 

S16 

Standard deviation of altitude, velocity, path 
angle, azimuth, latitude and longitude. Array of 
6 values. Input when lERROR « 3. 

SI 

Nominal state parameters at last stage burnout If 
covariance matrix propagation is not selected. 

If propagation Is selected, SI Is state at 
burnout of next to last stage. Array of 6 
values. Order Is altitude, velocity, path angle, 
azimuth, latltxide and longitude. 

S2 

Nominal state parameters at last stage burnout. 
Input If TCOAST Is non-zero. Array of 6 values. 
Order Is altitude, velocity, path angle, azimuth, 
latitude and longitude. 
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NAMELIST Input Definitions (Continued) 


TCOAST Nominal coast time to propagate covariance 

matrix. Input zero if propagation is not 
desired. (0. built-in) 

TITLEl Title information printed at top of each page of 

TITLE2 output. 72 characters maximum each. 

3.4 Output Description 

Both the fixed field and NAMELIST input data are listed verbatim on the 
output listing. These lists provide a quick check of the input data for 
format correctness and validity. 

Following the input data lists, the flight experience data are provided in 
a labeled format. The next page provides the flight experience data after the 
altitude, velocity, path angle errors have been normalized and the latitude, 
longitude errors have been adjusted. 

The following page provides error statistics for only those flights with 
latitude and longitude errors. The next page provides error statistics based 
upon all flight samples provided. The altitude, velocity, path angle and 
azimuth error statistics are obtained from those flight samples identified for 
that purpose. The latitude and longitude error statistics are derived from 
only those flights available for latitude and longitude statistics. Unless 
the covariance matrix is propagated, the final error covariance matrix is 
included on this page. Following the covariance matrix, the sensitivity 
matrix is provided, which when properly sampled will produce the covariance 
matrix. 

If the covariance matrix is propagated, additional matrices are shown on 
the next two pages. The sensitivity matrix obtained from the final error 
covariance matrix is also included. 
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TABLE 3.1 


SCOUT FLIGHT EXPERIENCE 
AT ORQIT INSERTION 


ERROR FRON PREDICTED 


VEHICLE 

INSERTION 



PATH 





ALTITUDE 

ALTITUDE 

VELOCITY 

ANGLE 

AZIMUTH 

LATITUDE 

LONGITUDE 


N.NI. 

N.HI. 

FPS 

DEG 

DEG 

DEG 

DEG 

S-136 

562.6 

9.7 

-58.5 

.77 

.28 



S-131 

612.5 

6.9 

-280.0 

— .36 

-.15 



S-138 

383.8 

3.2 

-107.9 

.15 

.77 



S-139 

40A.3 

8.9 

-29.2 

.05 

.19 



S-UO 

506.8 

-.2 

67.2 

-.37 

-.91 



S-1A2 

487.1 

4.1 

142.9 

— .86 

-.29 



S-1A3 

462.2 

4.2 

112.3 

.14 

-.28 



S-1A5 

200.1 

-3.5 

15.0 

-.75 

.48 



S-146 

495.2 

-3.7 

-24.3 

-.63 

0.00 



S-1A7 

350.6 

3.8 

11.1 

-.98 

-.24 



S-1A8 

197.7 

2.8 

-31.4 

-.17 

-.55 



S-1A9 

576.3 

0.0 

-38.6 

.16 

-1.18 



S-150 

178.5 

1.0- 

-31.1 

-.63 

-.14 

.0638 

.0339 

• S-154 

577.0 

-.7 

-57.5 

.12 

.26 

-.0473 

-.0361 

S-153 

117.1 

-.5 

-56.9 

-.63 

-.13 



S-155 

279.4 

-5.5 

94.3 

-.06 

.18 

-.1042 

.0115 

S-156 

592.7 

2.0 

-53.4 

-.05 

-.43 

-.0108 

.3134 

S-157 

565.9 

2.8 

-34.4 

-.15 

-.73 

-.0429 

.3543 

S-158 

235.3 

-2.1 

8.0 

-.17 

.70 

-.0729 

-.0749 

S-162 

581.2 

-1.7 

-29.1 

-.41 

-.07 

-.2190 

.4662 

S-161 

190.5 

-5.8 

17.0 

-.46 

-1.08 

-.0774 

.0137 

S-165 

376.8 

-2.3 

72.4 

-.72 

-1.36 

-.0345 

.1321 

S-167 

147.2 

.4 

30.7 

-.57 

-.25 

-.0300 

.0690 

S-172 

216.0 

-5.2 

-91.1 

-.13 

-.93 

-.0601 

.0344 

S-169 

214.7 

-1.3 

-35.0 

-.13 

.31 

-.0152 

.0477 

S-176 

588.3 

-10.1 

2.0 

-1.01 

.02 

-.0233 

.1267 

S-17A 

169.9 

-3.1 

-36.2 

-.17 

-.69 

-.0656 

-.0159 

S-175 

294.4 

-.7 

7.6 

-.07 

.11 

-.2135 

.0215 

S-173 

115.5 

4.5 

-93.0 

-.08 

-1.05 



S-177 

323.3 

• 6 

-79.4 

-.63 

-.54 

.2250 

.0386 

S-180 

486.0 

4.7 

-202.0 

-.24 

.06 

.2448 

-.1021 

S-163 

120.0 

-.4 

-129.0 

.49 

-.12 

.1882 

.7061 

S-163 

297.0 

-9.2 

47.0 

-.47 

0.00 

.0015 

.0170 

S-18A 

263.5 

6.7 

-8.0 

.26 

.32 

-.1355 

-.0692 

S-182 

449.6 

8.0 

-177.0 

-.02 

.13 

.2266 

.1281 

S-170 

299.7 

0.0 

-29.5 

-.76 

-.60 

.2533 

-.0426 

S-185 

151.2 

-3.5 

70.5 

-1.14 

.34 

-.0743 

-.0523 

S-181 

129.5 

,-3.4 

84.2 

-.41 

-.28 

-.0777 

.0682 

S-178 

604.4 

14.5 

-220.2 

-.03 

.17 

.3180 

.0550 

S-190 

123.1 

2.4 

58.2 

.03 

.07 

-.0017 

.0943 

S-188 

405.0 

-2.9 

149.5 

-.47 

-.62 

-.1845 

.1735 

S-101 

215.2 

.6 

45.7 

-.56 

-.20 

-.0857 

-.2038 

S-186 

124.2 

-4.2 

3.9 

-.24 

.69 

-.0147 

-.5603 

S-189 

213.5 

-5.0 

61.1 

-.05 

.28 

.0101 

-.3837 

S-197 

270.2 

4.0 

39.1 

.15 

0.00 

.0601 

.1696 

S-19A 

271.0 

.8 

14.1 

.03' 

-.43 

.1489 

.0124 

S-1R5 

194.8 

3.4 

-25.8 

-.09 

.77 

.0623 

-.1025 


NOTEl LATITUDE# LONGITUDE ERRORS NOT SHORN ARE NOT AVAILABLE 
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4.0 

This sectioa provides a brief description of each subroutine of STEP. 

4.1 STEP (Main Program) 

The main program initializes the input data defaults; calls the two input 
subroutines; normalizes the errors in altitude, velocity and flight path 
angle; adjusts the errors in latitude and longitude; calculates the error 
covariance matrix at the input epoch; calculates the sensitivity matrix; 
calculates the boost covariance matrix; propagates the covariance matrix; and 
outputs the results. 


4.2 CAROS 

Subroutine CAROS reads the fixed field data in alphanumeric format; writes 
the data as read on the output file in alphanumeric format; and writes the 
data as read on Unit 8 for subsequent reading by the main program in floating 
point format. CARDS counts the number of samples read and writes error 
messages if there are no samples or if they exceed the maximum of 100. 


4.3 CONIC 

Subroutine CONIC initializes a conic path from an input trajectory state 
(altitude, velocity, flight path angle, azimuth, latitude and longitude) for 
subsequent propagation of the state along the conic path by an input time 
increment. CONIC verifies that the conic path is elliptical and, if so, 
calculates Che orbital elements and coordinate transformations from the 
spherical state Co inertial cartesian components . CONIC is called one time 
per conic path. 
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4.4 CORCO 


Subroutine CORCO calculates mean values, standard deviations and the 
correlation coefficient of two independent variables from random samples of 
each. On option, the mean values may be set to zero. 


4.5 COVR 

Subroutine COVR generates a symmetric error covariance matrix of six 
parameters. COVR is called as each set of error samples is generated. COVRl 
(an entry point) is called after all sets of samples have been generated and a 
covariance matrix is desired. 


4.6 EIGEN 

Subroutine EIGEN calculates the eigenvalues and eigenvectors of a real 
symmetric matrix. These values are used In STEP to generate a sensitivity 
matrix of six pseudo. Independent, uncorrelated error sources to the six 

m 

trajectory parameters. 


4.7 INIT 

Subroutine INIT Initializes constants used by several of the subroutines, 
which are available to the subroutines via labeled common DIG. Entry point 
INITl Initializes the parameters used to obtain random numbers. These 
parameters are available to the subroutines via labeled common BLK4. 

4.8 INPUT 

Subroutine INPUT reads Input data cards In a modified NAMELIST format. 
Titling Information on title cards Is placed In appropriate arrays for use 
by the main program. Non-title cards are written on unit 8 for a NAMELIST read 
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by the main program 


4.9 INTER 

Subroutine INTER is a second-order interpolater of two variables. It 
selects the four closest data points to the desired value of the independent 
variable and interpolates or extrapolates for the value of the dependent 
variable. 


4.10 NEWTON 

Subroutine NEWTON iterates for the eccentric anomaly corresponding to a 
value of time along a conic. If the interation fails, a diagnostic is 
written. If the iteration is successful, the radius and true anomaly are 
calculated. 


4.11 NORRAN 

Subroutine NORRAN generates a normally distributed random number from the 
set of numbers which have a mean value and standard deviation as supplied to 
the subroutine. 


4.12 TSTEP 

Subroutine TSTEP propagates a trajectory state along a conic by a given 
time increment. Subroutine CONIC is used to initialize the conic from 
the initial state parameters. Subroutine NEWTON is used to iterate on 
eccentric anomaly. TSTEP updates- the trajectory state parameters at the new 
time. TSTEP is called each time the trajectory state is to be propagated. 
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4.13 DNIRAN 


Subroutine UNIRAN generates a unlfonoally distributed random number between 
zero and one. This remdom number is used when analyzing tipoff error sources* 
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5.0 PROGRAM CODIHG 


This section presents details about the progreun coding. Included are 
flowcharts of each subroutine, FORTRAN listings of each subroutine and 
definitions of the FORTRAN variables. The information presented in this 
section is intended to be helpful in developing a thorough understemding of 
STTEP and in making modifications to the program. 

5.1 Subroutine Flowcharts 

Flowcharts are presented in Figures 5.1 through 5.9. The flowchart 
conventions used are defined in Section 2.0 of this report. 

5.2 FORTRAN Listings 

STEP is coded in FORTRAN IV, Reference (2), on the GDC GIBER 175 computer 
with the NOS/BE 1.4 operating system. Listings of the FORTRAN coding are 
presented in Appendix B. 


5.3 FORTRAN Variable Definition 

Definitions of the FORTRAN variables are presented below. This information 
is normally used only when making modifications to the program. 

ALT Ncmiinal altitude, n.mi. 

BEST Array equivalenced to SENl, SEN2, SEN3, SEN4 

CE East component of position error, ft 

GLAT Gosine .of nominal latitude 

GMATl Error covariance matrix of last stage boost 

CMAT2 Error covariance matrix eifter propagation 
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CN 

CON 

COV 

CPI 

CP2 

CP3 

CZ 

C2ET 

D 

DATAG 

DATAH 

DATAV 

DCR 

DELT 


DEVGl 

DEVG2 

DEVHl 

DEVH2 

DEWl 

DEW2 

DGl 

DG2 

DHl 

OH2 

DLGl 

DLL 

DR 

DVl 

DV2 

DX 

DZl 

DZ2 

EHAG 

EPAR 


North component of position error, ft 
Radians per degree, .0174532925 
Input data 

Working arrays for correlation coefficients 

Cosine of azimuth error 
Cosine of nominal azimuth 
Working array 

Input data for normalizing Input flight errors 

Crossrange error, n.ml. 

Time Increment for propagation of state 
parameters, sec 

Interpolated values of DATAG 

Interpolated values of DATAH 

Interpolated values of DATAV 

Path angle errors of flight sample 

Altitude errors of flight sample 

Longitude error of flight sample 
Latitude error of flight sample 
Range error, n.ml. 

Velocity errors of flight sample 

Working array of trajectory parameter errors 
Azimuth errors of flight sample 

Input data 

Array of normally distributed random numbers 
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EXIC 

FT 

GAUSS 

GN 

HINJ 

ICNT 

ICODE 

lERHOR 

IWOHO 

KERR 

MEAN 

MEANG 

HEANH 

MEANL 

MEANLG 

MEMV 

MEANZ 

N 

NDATA 

NERROR 

NFliAG 

NORM 

NFROB 

NSAMP 

N1 

N2 


Working array of trajectory pareuneter errors 

Feet per n.mi. , 6076.11549 

Normally distributed random number 

Earth gravitational constant, 1.4076576E16 
f t3/sec2 

Random numbers used in random number generator 
Input data 

Counter used in sxibroutine NORRAN 
Value of flag input with flight sample 
Input data 

Array of flags input with flight seuaples 
Diagnostic flag set in subroutine NEWTON 
Flag for calculating mean values 
Mean value of path angle, deg 
Mean value of altitude, n.mi. 

Mean value of latitude, deg 
Mean value of longitude, deg 
Mean value of velocity, fps 
Mean value of azimuth, deg 

Number of flight samples 

Number of data values of HINJ 

Input data 

Working flag 

Input data 

Case number 

Input data 

Number of flight samples with 
ICODE - 1 or 2 

Number of flight samples with 
ICODE - 2 or 3 
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N3 

PI 

R 

RAD 

RE 

RHO 

RHOGL 

RBOGLG 

RHOGZ 

RHCHG 

RHOHI> 

RHOHLG 

RHOHV 

RHOHZ 

RH0LL6 

RHOVG 

RHOVL 

RHOVLG 

RHOVZ 

RHOZL 

RHOZLG 

ROTATE 

RR 

SC 

sEta 

SEN2 

SER3 

SEN4 


Rinnber of flight saonples with 
ICODE = 1, 2 or 3 

3-1415927 

Eigenvectors 

Degrees per radian, 57.2957795 
Earth equatorial radius, 20925741 ft 
Input data 

Correlation coefficient of path angle and latitude 

Correlation coefficient of path angle and 
longitude 

Correlation coefficient of path angle and azimuth 
Correlation coefficient of altitude and path angle 
Correlation coefficient of altitude and latitude 
Correlation coefficient of altitude and longitude 
Correlation coefficient of altitude and velocity 
Correlation coefficient of altitude and azimuth 
Correlation coefficient of latitude and longitude 
Correlation coefficient of velocity and path angle 
Correlation coefficient of velocity and latitude 
Correlation coefficient of velocity and longltvule 
Correlation coefficient of velocity and azimuth 
Correlation coefficient of azimuth and latitude 
Correlation coefficient of azimuth and longitude 
Earth rotation rate, 7.29211B-5 radians per sec. 
Eqpilvalenced. to R 
Sensitivity coefficients 
Input data 
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SI 6 

SIGALT 

SIGG 

SIGH 

SIGL 

SIGLG 

SIGM& 

SIG7 

SIGZ 

SS 

STATE 

SUMSQl 

ST3MSQ2 

SUMSQ3 

SUMl 

SXM2 

SUM3 

SZ 

SZET 

51 

52 
T 

TCALL 

TCOAST 

THBT 

TITUEl 

TITLE2 

TO MAG 

TZERO 


Input data 

Standard deviation of altitude, n.mi. 

Standard deviation of path angle, deg 

Standard deviation of altitude, ft 

Standard deviation of latitude, deg 

Standard deviation of longitude, deg 

Standard deviation of normally distributed 
random numbers* Used in subroutine NORHAM 

Standard deviation of velocity, fps 

Standard deviation of azimuth, deg 

Input state SI propagated b7 TCOAST 

Trajectory state to be propagated 

Working arrays of stmnnatlon of errors squeured 

Working arrays of summation of errors 

Sine of azimuth error 
Sine of nominal azimuth 
Input data 
Input data 

Time of propagated state referenced to TZERO, sec* 

Time of propagated state referenced to perigee, 
sec* 

Input data 

Random direction of tipoff, deg 
Input data 

Random value of tipoff magnitude, deg. 

Initial time of state to be propagated, sec. 
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w 

WORD 

X 

XAIiT 

XG 

XH 

XL 

XIAT 

XLG 

XHD 

XV 

XZ 

ZET 

ZRAN 


Wor]clng array for sensitivity matrix 
Array of flight sample data 
Random values of error sources 
Nominal altitude of flight sample, n.mi« 

Path angle error of flight seui^le, deg. 

Altitude error of flight sample, ft 
Latitude error of flight sample, deg. 

Nominal latitude of flight sample, deg. 

Longitude error of flight sample, deg. 

Mean value of normally distributed rzmdom ntimbers 
Velocity error of flight sample, fps 
Azimuth error of flight sample, deg. 

Nominal azimuth, deg. 

Uniformly distributed random number 


( 



Figure 5.1 

FLOWCHART OF MAIN PROGRAM S 












Figure 5.2 

FLOWCHART OF SUBROUTINE CARDS 





Figure 5.3 

FLOWCHART OF SUBROUTINE CONIC 
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Figure 5.4 

aOWCHART OF SUBROUTINE CORCO 
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Figure 5.5 

FLOWCHART OF SUBROUTINl 


( Enter 
COVR 


> 


Accumulate 
summations of . 
errors and 
errors squared 


( Enter 
COVRI 



Calculate 

standard deviations 
and correlation 
coefficients 
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Calculate 
covari ance 
matrix 



Return 








Figure 5.6 

FLOWCHART OF SUBROUTINE INPUT 







Figure 5.7 

FLOWCHART OF SUBROUTINE INTER 



Interpolate for 
dependent variable 
according to Index 
selected 


Select interpolator 
index to the 4 
closest values to 
the desired 
Independent variable 


Return 




F-igure 5.8 

FLOWCHART OF SUBROUTINE NEWTON 
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Figure 5.9 

FLOWCHART OF SUBROUTINE TSTEP 







HSFERENCES 


1. Vought Corporation Report No. 23.478/ "The Statistical Orbital Analysis 
Routine" dated 19 March 1971. 

2. Control Data Corporation/ "FORTRAN Extended Version 4 Reference 
Manual," Revision C, dated 15 April 1977. 


- 37 - 



APPENDIX A 
SAMPLE DATA CASES 


46 FLIGHT SAMPLES AT STAGE 4 BURNOUT 


S-136 

3. 

562.8 

9.7 

-58.9 

0.77 

0.28 




S-131 

4. 

612.5 

6.9 

-280.0 

-0.36 

-0.19 




S-138 

3. 

383.8 

3.2 

-107.9 

0.15 

0.77 




S-139 

3. 

404.3 

8.9 

-29.2 

0.05 

0.19 




S-140 

3. 

506.8 

-0.2 

67.2 

-0.37 

-0.91 




S-142 

3. 

487.1 

4.1 

142.9 

-0.86 

-0.29 




S-143 

3. 

482.2 

4.2 

112.3 

0.14 

-0.28 




S-145 

3. 

200.1 

-3.5 

15.0 

-0.75 

0.48 




S>146 

3. 

495.2 

-3.7 

-24.3 

-0.63 

0.0 




S-147 

3. 

390.6 

3.8 

11.1 

-0.98 

-0.24 




S-14B 

3. 

197.7 

2.8 

-31.4 

-C.17 

-0.59 




S-149 

3. 

576.3 

0.0 

-38.6 

0.16 

-1.18 




S-190 

2. 

178.9 

1.0 

-31.1 

-0.63 

-0.14 0.0638 

0.0339 

20.7161 

171.556 

S-154 

2. 

577.0 

-0.7 

-57.9 

0.12 

0.26 -0.0473 

-0.0361 

12.1787 

180.000 

S-153 

3. 

117.1 

—0.5 

-56.9 

-0.63 

-0.13 




S-155 

2. 

279.4 

-5.9 

94.3 

-0.06 

0.18 -0.1042 

0.0119 

19.0829 

169.640 

S-156 

2. 

582.7 

2.0 

-93.4 

-0.09 

-0.43 -0.0108 

0.3134 

12.0870 

180.000 

S-157 

2. 

565.9 

2.8 

-34.4 

-0.15 

-0.73 -0.0429 

0.3943 

12.1900 

180.001 

S-198 

2. 

235.3 

-2.1 

8.0 

-C.17 

0.70 -0.0729 

-0.0749 

18.2632 - 

179.999 

S-162 

2. 

581.2 

-1.7 

-29.1 

-0.41 

-0.07 -0.2190 

0.4662 

12.2504 

180.000 

S-161 

2. 

190.5 

-9.8 

17.0 

-0.46 

-1.08 -0.0774 

0.0137 

22.0366 

188.853 

S-165 

2. 

376.8 

-2.3 

72.4 

-0.72 

-1.36 -0.0349 

0.1321 

13.4573 

171.761 

S-167 

2. 

147.2 

0.4 

30.7 

-0.97 

-0.25 -0.0300 

0.0690 

23.1376 

184.350 

S-172 

2. 

216.0 

-9.2 

-91.1 

-0.13 

-0.93 -0.0801 

0.0344 

20.8663 

175.719 

S-169 

2. 

214.7 

-1.3 

-39.0 

-C.13 

0.31 -0.0152 

0.0477 

17.7236 

193.313 

S-176 

2* 

588.3 

-10.1 

2.0 

-1.01 

0.02 -0.0233 

0.1267 

11.9576 

180.000 

S-174 

2. 

169.9 

-3.1 

-38.2 

-0.17 

-0.69 -0.0656 

-0.0159 

39.9464 

103.438 

S-175 

2. 

294.4 

-0.7 

7.6 

-0.07 

0.11 -0.2135 

0.0215 

-2.7201 

88.955 

S-173 

3. 

115.5 

4.9 

-93.0 

-0.08 

-1.09 




S-177 

2. 

323.3 

0.6 

-79.4 

—0.63 

-0.54 0.2250 

0.0386 

23.3147 

137.246 

S-180 

2. 

486.0 

4.7 

-202. 

-0.24 

0.06 0.2448 

-0.1021 

23.9077 

135.496 

S-163 

2« 

120.0 

-0.4 

-129. 

0.49 

-0.12 0.1882 

0.7061 

-3.3437 

90.869 

S-18 3 

2. 

297.0 

-9.2 

47.0 

-0.47 

0. 0.0019 

0.0170 

15.4234 

172.737 

S-18 4 

2. 

263.9 

6.7 

-8.0 

0.26 

0.32 -0.1359 

-0.0692 

34.3083 

106.662 

S-1S2 

2. 

449.6 

8.0 

-177.0 

-0.02 

0.13 0.2266 

0.1281 

13.4634 

180.000 

S-170 

2. 

299.7 

-0.0 

-29.9 

-0.76 

-0.60 0.2933 

-0.0426 

-1.6043 

89.260 

S-18S 

2. 

191.2 

-3.5 

70.9 

-1.14 

0.34 -0.0743 

-0.0923 

22.3683 

180.860 

S-lBl 

2. 

129^.9 

-3.4 

84.2 

-0.41 

-0.28 -0.0777 

0.0682 

22.1681 

187.780 

S-178 

2. 

604.4 

14.9 

-220.2 

-0.03 

0.17 0.3180 

0.0950 

11.5952 

180.020 

S-190 

2. 

123.1 

2.4 

98.2 

0.03 

0.07 -0.0017 

0.0943 

-2.8471 

89.370 

S-188 

2. 

409.0 

-2.9 

149.5 

-0.47 

-0.62 -0.1849 

0.1739 

13.8902 

188.650 

S-191 

2* 

219.2 

0.6 

45.7 

-0.56 

-0.20 -0.0857 

-0.2038 

14.2879 

179.896 

S-186 

2. 

124.2 

-4.2 

3.9 

-0.24 

0.69 -0.0147 

-0.9603 

21.9731 

187.336 

S-1B9 

2. 

210.9 

-9.0 

61.1 

-0.09 

0.28 0.0101 

-0.3837 

13.0322 

187.980 

S-18 7 

2. 

270.2 

4.0 

39.1 

0.19 

0.00 0.0601 

0.1696 

-2.7540 

89.026 

S-194 

2. 

271.0 

0.8 

14.1 

C.03 

-0.43 0.1489 

0.0124 

-2.6965 

88.795 

S-195 

2. 

194.8 

3.4 

-25.8 

-0.09 

0.77 0.0623 

-0.1029 

19.0391 

179.991 


END 
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ROUTINE STEP 
DTD NAR 1981 
PROS. NO. 1 

TRAJECTORY TO 200 NH ALTITUDE 
WITH PROPAGATION 500 SAMPLES 


FLIGHTS FOR SCOUT FLIGHT EXPERIENCE ACCURACY 
standard DEVIATIONS 

ALT-VEL-GAM-IET-LAT-LONG (FT>FPS*DEG) 

« 

29297. 80.4 *402 *525 *1641 .2239 

4.822 NN 

MEAN VALUES 

- ALT-VEL-GAN-ZET-LAT-LONG 

2391. -11.5 -.260 >.152 >.0347 .0427 


CORRELATION COEFFICIENTS 
ALT>VEL-GAH>Z ET>LAT>LONG 

>.43278 .43300 .11934 .40168 .14018 

>.32800 >.10826 >.33822 >.20859 

.16758 >.21837 >.11681 

.07191 -.47921 

>.07939 


COVARIANCE MATRIX (A MATRIX) 

ALT VEL GAN ZET LAT LONG 

.85830E-»09 

>.10200E^07 .64719E>04 

• 51043E404 >.10618E4>02 .16191E-»^00 

.18355E+04 -.45721E+01 .35398E-01 .27558E+00 

.19316E4>04 -.44661E401 -.14423E-01 .61961E-02 .26943E>01 

.91966E>03 -.37578E401 -.10525E-01 >.56335E>01 -.29183E-02 .50149E-01 


EIGENVALUES 

• 85830E409 .52597E-»04 .28973E«00 .13160E-»00 .35253E-01 .11113E-01 

EIGENVECTORS 

.lOOOOE^Ol .11884E-02 >.21914E>05 >.39473E>05 >.11341E>06 -.31095E-05 

>.11884E>02 .10000E401 .44569E-03 .61192E-03 .39731E-03 .80723E-03 

.59470E-05 >.86539E>03 .16123E«00 .95227E«00 .40492E-02 .23910E>00 

• 21385E-05 >.4S4S8E>03 .95546E«00 -.17894E-»00 .22700EO0 .59516E-01 

.22503E-05 >.41268E>03 >.72S90E>02 -.23903E-^00 >.39090E4^00 .88883E400 

.10715E-05 >. 5066SE-03 >.24707E-»00 -.63533E-01 .89199E^00 .37319E>00 

SENSITIVITY COEFFICIENTS 
ALT>VEL>GAH-ZET>LAT>LONG (TOP TO BOTTOM) 

.29297E+05 .86188E-01 -.11795E-05 -.14319E-05 -.21294E-07 >.32780E>06 

>.34816E«02 .72524E+02 .23990E-03 .22198E-03 .74598E-04 .85096E-04 

.17423E«00 >.62761E>01 .8678SE-01 .34545EfOO .76026E-03 .27322E-01 

.62651E-01 -.32968E-01 .51429E*00 -.64912E-01 .4262ZEH)1 .62740E-02 

.65933E-01 -.29929E-01 >.39|073E>02 >.86710E>01 -.73394E-01 .93698E-01 

•31391E-01 -.36744E-01 >.13299E«00 -.23047E-01 .16748E+00 .39340E-01 
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I routine step 

OTO MAR 1961 
PROS* NO. 1 

TRAJECTORY TO 200 NN ALTITUDE 
WITH PROPAGATION 500 SAMPLES 


COVARIANCE MATRIX OF LAST STAGE BOOST 


ALT 

VEL 

GAM • ' 

ZET 

LAT 

LONG 

•4AA72E408 

-.97678E>05 

.41439E*03 





.610A6E-»04 

-.13451E402 

.83797E+00 




-.17395E+03 

63128E«01 

-.2447AE-01 

.76739E«00 



•15782E^02 

-.117A0E-01 

•21447E-02 

.44487E-02 

.57817E-04 


•.33073E-K)2 

. 30522 E-01 

•44994E-02 

•21484E-02 

•77457E-04 

•13641E-03 


SAMPLED COVARIANCE MATRIX 

AT INPUT 

EPOCH 



ALT 

VEL 

GAN 

ZET 

LAT 

LONG 

•91102E409 






-•13486E407 

« 76936 E«04 





• 55736E-»04 

-.13026E402 

• 17963E-«>00 




•19986E404 

-.12525E+01 

.25690E-01 

.27693E-^00 



• 24830E-»04 

-.68055E+01 - 

.14711E-01 

.97317E-02 

.31043E-01 


•98971E«03 

-.57254E+01 - 

.30337E-02 

-.61651E-01 

-.25342E-02 

•52863EH)1 


COVARIANCE 

MATRIX AFTER 

PROPAGATION 




ALT 

VEL 

GAN 

ZET 

LAT 

LONG 

.31132E«10 






-.44132E^07 

.12380E40S 





.18540E4-05 

-.20220E4-02 

.15571E«00 




.38874E>04 

-.43581E*01 

.22484E-01 

.22453E«00 



.47692E+04 

-.1S010E4^02 

-.57539E-02 

.12444E-01 

.50500E-01 


-.12004E4-04 

-.29157E+01 

-.20754E-01 

-.11936E«00 

-.18270E-02 

•10248E«00 
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' ROUTINE STEP 

OTD KAR 1981 
PROS. NO. 2 

TRAJECTORY TO 200 NH ALTITUDE 
NO PROPAGATION 500 SAMPLES 


NORNALIZEO DEVIATIONS AND ADJUSTED LATITUDE* LONGITUDE DEVIATIONS 




ALTITUDE 

• 90.0 

N Ml 






AZIMUTH- 

178.900 

DEG 






LATITUDE 

• 30.6500 

DEG 




SAMPLE 

ALT 

VEL 

GAMMA 

ZETA 

LAT 

LONG 

CODE 

NUMBER 

FT 

FPS 

DEG 

DEG 

DEG 

DEG 


• 1 

58938. 

-58.5 

.770 

.280 

0.0000 

0.0000 

3 

2 

19444. 

-107.9 

.150 

.770 

o.oooe 

0.0000 

3 

3 

54077. 

-29.2 

.050 

• 190 

0.0000 

0.0000 

3 

4 

-1215. 

67.2 

-.370 

-.910 

0.0000 

0.0000 

3 

3 

24912. 

142.9 

-.860 

-.290 

0.0000 

0.0000 

3 

6 

25520. 

112.3 

.140 

-.280 

0.0000 

O.OQOO 

3 

7 

-21266. 

15.0 

-.750 

.480 

0.0000 

0.0000 

3 

8 

-22482. 

-24.3 

-.630 

0.000 

0.0000 

0.0000 

3 

9 

23C89. 

11.1 

-.980 

-.240 

0.0000 

0.0000 

3 

10 

17013. 

-31.4 

-.170 

-.550 

0.0000 

0.0000 

3 

11 

0. 

-38.6 

.160 

-1.180 

0.0000 

0.0000 

3 

12 

6076. 

-31.1 

-.630 

-.140 

.0592 

.0460 

2 

• 13 

-4253. 

-57.5 

.120 

.260 

-.0480 

-.0400 

2 

14 

-3038. 

-56.9 

-.630 

-.130 

0.0000 

0.0000 

3 

15 

-33419. 

94.3 

-.060 

.180 

-.1046 

-.0066 

2 

16 

12152. 

-53.4 

-.050 

-.430 

-.0049 

.3564 

2 

17 

17013. 

-34.4 

-.150 

-.730 

-.0362 

.4034 

2 

18 

-12760. 

8.0 

-.170 

.700 

-.0743 

-.0810 

2 

19 

-10329. 

-29.1 

-.410 

-.070 

-.2102 

• 5344 

2 

20 

-35241. 

17.0 

-.460 

-1.080 

-.0740 

.0301 

2 

21 

-13975. 

72.4 

-.720 

-1.360 

-.0502 

.1432 

2 

22 

2430. 

30.7 

-.570 

-.250 

-.0238 

.0767 

2 

23 

-31596. 

-91.1 

-.130 

-.930 

-.0818 

.0321 

2 

24 

-7899. 

-35.0 

-.130 

.310 

-.0034 

.0555 

2 

25 

-61369. 

2.0 

-1.010 

.020 

-.0209 

.1446 

2 

26 

-18836. 

-38.2 

-.170 

-.690 

-.0039 

-.0776 

2 

27 

-4253. 

7.6 

-.070 

.110 

-.0217 

-.2481 

2 

28 

27343. 

-93.0 

-.080 

-1.050 

0.0000 

0.0000 

3 

29 

3646. 

-79.4 

-.630 

-.540 

.1446 

.2046 

2 

30 

28558. 

-202.0 

-.240 

.060 

.2422 

.1165 

2 

31 

-2430. 

-129.0 

• 490 

-.120 

-.6980 

• 2467 

2 

32 

-*5900. 

47.0 

-.470 

0.000 

-.0003 

.0191 

2 

33 

40710. 

-8.0 

.260 

.320 

.0131 

-.1703 

2 

34 

48609. 

-177.0 

-.020 

.130 

.2289 

.1397 

2 

35 

0. 

-29.5 

-.760 

-.600 

.0442 

.2941 

2 

36 

-21266. 

70.5 

-1.140 

.340 

-.0759 

-.0532 

2 

37 

-20659. 

84.2 

-.410 

-.260 

-.0670 

.0865 

2 

38 

88104. 

-220.2 

-.030 

.170 

.3190 

.0554 

2 

39 

14583. 

58.2 

.030 

.070 

-.0942 

-.0011 

2 

40 

-17621. 

149.5 

-.470 

-.620 

-.1533 

.2293 

2 

41 

3646. 

45.7 

-.560 

-.200 

-.0891 

-.2278 

2 

42 

-25520. 

3.9 

-.240 

.690 

-.0908 

-.5949 

2 

43 

-30381. 

61.1 

--.050 

.280 

-.0490 

-.4309 

2 

44 

24304. 

39.1 

.150 

0.000 

-.1693 

.0703 

2 

45 

4861. 

14.1 

.030 

-.430 

-.0127 

.1731 

2 

46 

20659. 

-25.8 

-.090 

.770 

.0604 

-.1164 

2 



ROUTINE STEP 
DTO NAR 1981 
PROS. NO. 2 


TRAJECTORY TO 200 NN ALTITUDE 
NO PROPAGATION 300 SAMPLES 


FLIGHTS WITH LATITIUOE/LONGITUDE DEVIATIONS 
33 FLIGHT SAMPLES 

STANDARD DEVIATIONS 

ALT-VEL-GAM-ZET-LAT-LONG (FT»FPS.DEGI 


29603. 83. A .361 .503 .16A1 .2239 

4.872 NM 


MEAN VALUES , 

ALT-VEL-GAM-ZET-LAT-LONG 

-2799. -13.2 -.263 -.123 -.0347 .0427 


CORRELATION COEFFICIENTS 
ALT-VEL-GAM-ZET-LAT-LONG 


-.56673 


.40712 

-.30939 


.16499 

-.08431 

.25407 


.40168 

-.33822 

-.21837 

.07191 


.14018 

-.20859 

-.11681 

-.47921 

-.07939 


ALT 


COVARIANCE MATRIX (A MATRIX) 

VEL GAM ZET LAT 


LONG 


.87636E409 

.13999E+07 

.43566E^Q4 

.24649E-<’04 

.19318E>04 

.92929E>03 


.69624E«04 

-.93318E-»01 

-.33585E-»01 

-.46323E+01 

-.38976E4^01 


•13067E400 

.46347E-01 

-.12957E-01 

-.94536E-02 


.25467E-f00 

.39564E-02 

-.54136E-01 


.26943E-01 
• -.29183E-02 


.50149E-01 


-A5- 



TRAJECTORY TO 200 NH ALTITUDE 
WITH PROPAGATION 900 SAMPLES 


FLIGHTS WITH LATITIUDE/LOMGITUDE DEVIATIONS 
33 FLIGHT SAMPLES 

STANDARD DEVIATIONS 

ALT-VEL-GAM-ZET-LAT-LONG (FT«FPS«DEG) 


29603. 83.4 ,361 

4.872 NM' 

.909 

.1641 

.2239 

MEAN VALUES 

alt-vel-gam-zet-lat-long 

-2799, -13.2 -.265 

-.123 

-.0347 

• 0427 

CORRELATION COEFFICIENTS 
ALT-VEL-GAM-ZET-LAT-LONG 




-.56673 .40712 .16499 

-.30939 -.08451 

•29407 

.40168 

-.33822 

-.21837 

.07191 

.14018 

-.20859 

-.11681 

-,47921 

-.07939 



ALT 


COVARIANCE MATRIX (A MATRIX) 

VEL GAM ZET LAT . 


.87636E«09 
.13999E+07 
• 43566E4^0A 
.24649E4^04 
•195I8E«04 
.92929E403 


• 69624E-»04 
93318 E+01 
-.35505E+OI 
-.46323E+01 
-.38976E+01 


.X3C67E-»00 

•46347E>01 

-.12937E-01 

-.94956E-02 


.29467E>00 

•99964E~02 

-.54156E-01 


.26943E-01 

-.29183E-02 


ROUTINE STEP 
OTO MAR 1981 
PROB. NO. 1 


LONG 


•90149E-01 


ROUTINE STEP 
DTD MAR 1981 
PROS. NO. 1 


TRAJECTORY TO 200 NN ALTITUDE 
IfITH PROPAGATION 500 SAMPLES 


COVARIANCE MATRIX AT LAST STAGE BURNOUT 


ALT 

.31576E410 

.45109E+07 

.24645E«05 

.371146^4 

•47850E404 

•11674E^04 


VEL 

.12995E405 
-.33671E«02 
-.10671E*0Z 
>. 15021 E-»02 
>.28852E*01 


GAM 


«99368E«00 

-.19893E-02 

>.36091E-02 

-.16255E-01 


ZET 


.10119E401 

•16893E>01 

-.11721E+00 


lAT 


. 50558 E*01 
-.17495E-02 


EIGENVALUES 

.31576E-t-10 .65504E-»04 .10225E*01 .79928E-H)0 .83924E«01 

EIGENVECTORS 


• 1C000E-»01 
.14286E-02 
.78047E-05 
•11754E-05 
.15154E-05 
«36969E>06 


. 14286 E-02 
.10000E<^D1 
. 23448 E-03 
-.81974E-03 
-.12497E-02 
-.69504E-03 


.88576E-06 

.76165E-03 

-.13018E«00 

.98333E-»00 

•10391E-01 

-.12653E400 


-.81140E-05 

-.21010E-03 

•99012E400 

•12774E+00 

-•49461E-01 

-.30003E-01 


.12497E-05 

.65864E-03 

.eC103E«02 

.12904E400 

-.10453E400 

.98608E>00 


SENSITIVITY COEFFICIENTS 
ALT-VEL-6AH-ZET-LAT-L0NG (TOP TO BOTTOM) 


•56193E>05 

.e0276E«02 

.43857E*00 

.66048E-01 

.85154E-01 


. 11562 E^OO 
•e0934E«02 
. 18978 E-Ol 
-.66346E-01 
10114E-»00 


.89567E-06 -.72542E-05 .36203E»06 

•77017E-03 -.18784E-03 .19081E-03 

-.13163E«00 .88519E-»00 ' .23206E-02 

.99434E«00 .11420E-f00 .37384E-01 

.10S06E-01 -.44220E-01 -.3C283E-01 


.20774E-01 -.56253E-01 


-«12794E*00 -.26S23Er01 .28566E-»00 


LONG 


.10262E«00 


.30506E-01 


-.10137E-07 

.13090E-02 

.51511E-01 

.96536E-02 

.99324E400 

.10361E4^00 


-.17704E-08 
•22B64E-03 
.89969E-02 
.16861E-02 
• 17348E-»0C 
.18096EHU 



ROUTINE STEP 
OTO NAR 1981 
PROB. MO. 1 


TRAJECTORY TO 200 NN ALTITUDE 
MXTH PRQPA&ATION 300 SAMPLES 


NORMALIZED DEVIATIONS AND ADJUSTED LATITUDE* LONGITUDE DEVIATIONS 




ALTITUDE 

• 90. 0 

N MI 





■ AZIMUTH- 

178.900 

DEG 





LATITUDE 

• 30.6500 

DEG 



SAMPLE 

ALT 

VEL 

GAMMA 

ZETA 

LAT 

LONG 

NUMSER 

FT 

FPS 

DEG 

DEG 

DEG 

DEG 

1 

58938* 

-58.5 

.770 

.280 

0.0000 

0.0000 

2 

19999. 

-107.9 

.150 

.770 

0.0000 

0.0000 

3 

59077. 

-29.2 

.050 

.190 

0.0000 

0.0000 

A 

•1215. 

67.2 

-.370 

-.910 

0.0000 

0.0000 

5 

29912. 

192.9 

-.860 

-.290 

0.0000 

0.0000 

6 

25520. 

112.3 

.190 

-.280 

0.0000 

0.0000 

7 

-21266* 

15.0 

-.750 

*980 

0.0000 

0.0000 

8 

-22982. 

-29.3 

-.630 

0.000 

0.0000 

0.0000 

9 

23089. 

11.1 

-.980 

-.290 

0.0000 

0.0000 

10 

17013. 

-31.9 

-.170 

-.550 

0.0000 

0.0000 

11 

0. 

-38.6 

.160 

-1.180 • 

0.0000 

0.0000 

12 

6076. 

-31.1 

-.630 

-.190 

.0592 

• 0960 

13 

-9253. 

-57.5 

.120 

.260- 

-.0980 

-.0900 

19 

-3038. 

-56.9 

-.630 

-.130 

0.0000 

0.0000 

15 

-33919. 

99.3 

-.060 

.180 

-.1096 

-.0068 

16 

12152. 

-53.9 

-.050 

-.930 

-.0099 

.3569 

17 

17013. 

-39.9 

-.150 

-.730 

-.0362 

.903 9 

18 

-12760. 

8.0 

-.170 

.700 

-.0793 

-.0810 

19 

-10329. 

-29.1 

-.910 

-.070 

-.2102 

.5399 

20 

-35291. 

17.0 

-.960 

-1.080 

-.0790 

.0301 

21 

-13975. 

72.9 

-.720 

-1.360 

-.0502 

.1932 

22 

2930. 

30.7 

-.570 

-.250 

-.0238 

.0767 

23 

-31596. 

-91.1 

-.130 

-.930 

-.0818 

.0321 

29 

-7899. 

-35.0 

-.130 

• 310 

-.0039 

.0555 

25 

-61369. 

2.0 

-1.010 

.020 

-.0209 

.1996 

26 

-18836. 

-38.2 

-.170 

-.690 

-.0039 

-.0776 

27 

-9253. 

7.6 

-.070 

.110 

-.0217 

-.2981 

28 

27393. 

-93.0 

-.080 

-1.050 

0.0000 

o.cooo 

29 

3696. 

-79.9 

-.630 

-.590 

.1996 

.2096 

30 

28558. 

-202.0 

-.290 

.060 

.2922 

.1165 

31 

-2930. 

-129.0 

.990 

-.120 

-.6980 

.2967 

32 

-55900. 

97.0 

-.970 

0.000 

-.0003 

.0191 

33 

90710. 

-8.0 

.260 

• 320 

.0131 

-.1703 

39 

98609. 

-177.0 

-.020 

.130 

.2289 

.1397 

35 

0. 

-29.5 

-.760 

-.600 

.0992 

.2991 

36 

-21266. 

70.5 

-1.190 

.390 

-.0759 

-.0532 

37 

-20659. 

89.2 

-.910 

-.280 

-.0670 

.0865 

38 

88109. 

-220.2 

-.030 

• 170 

.3190 

.0559 

39 

19583. 

58.2 

.030 

.070 

--0992 

-.coil 

90 

-17621. 

199.5 

-.970 

-.620. 

-.1533 

.2293 

91 

3696. 

95.7 

-.560 

-.200 

-.0891 

-.2278 

92 

-25520. 

3.9 

-.290 

.690 

-.0908 

-.5999 

93 

-30381. 

61.1 

-.050 

.280 

-.0990 

-.9309 

99 

29309. 

39.1 

.150 

0.000 

-.1693 

.0703 

95 

9861. 

19.1 

.030 

-.930 

-.0127 

.1731 

96 

20659. 

-25.8 

-.090 

.770 

.0609 

-.1169 


CODE 


-A8- 


NNIMIVMNNMMIVMMIMIMIMMIUMOtMNMIVtNritNINtNINtNNIMWNIViWbIbIWWWWWWWUt 



ROUTINE STEP 
DTD MAR 1981 
PROD. NO. 1 

TRAJECTORY TO 200 NN ALTITUDE 
UlTH PROPAGATION SOO SAMPLES 


PLIGHT. RESULTS AT INPUT TINE 


SAMPLE 

VEHICLE 

PREDALT 

DELTA 

DELTA 

DELTA 

DELTA 

DELTA 

DELTA 

NUMBER 

NUMBER 

- 

ALT 

VEL 

GAM 

AZ 

LAT 

LONG 



N.MI. 

N.MI. 

’fps 

DEG 

OEG 

OEG 

OEG 

1 

S-136 

962.8 

9.700 

-58.9 

.770 

.280 

0.0000 

0.0000 

2 

S-131 

612.5 

6.900 

-280.0 

-.360 

-.190 

0.0000 

0.0000 

3 

S-138 

383.8 

3.200 

-107.9 

.150 

.770 

0.0000 

0.0000 

A 

S-139 

404.3 

8.900 

-29.2 

.050 

.190 

0.0000 

0.0000 

9 

S-140 

506.8 

-.200 

67.2 

-.370 

-.910 

0.0000 

0.0000 

6 

S-142 

487.1 

4.100 

142.9 

-.860 

-.290 

0.0000 

0.0000 

7 

S-143 

482.2 

4.200 

112.3 

• 140 • 

-.280 

0.0000 

0.0000 

8 

S-149 

200.1 

-3.500 

15.0 

-.750 

.480 

0.0000 

0.0000 

9 

S-146 

495.2 

-3.700 

-24.3 

-.630 

0.000 

0.0000 

0.0000 

10 

S-147 

350.6 

3.800 

11.1 

-.980 

-.240 

0.0000 

0.0000 

11 

S-148 

197.7 

2.800 

-31.4 

-.170 

-.590 

0.0000 

0.0000 

12 

S-149 

576.3 

0.000 

-38.6 

.160 

-1.180 

0.0000 

0.0000 

13 

S-190 

178.5 

1.000 

-31.1 

-.630 

-.140 

.0638 

.0339 

14 

S-194 

977.0 

-.700 

-57.9 

.120 

.260 

-.0473 

-.0361 

15 

S-193 

117.1 

-.500 

-56.9 

-.630 

-.130 

0.0000 

0.0000 

16 

S-199 

279.4 

-9.500 

94.3 

-.060 

.180 

-.1042 

.0115 

17 

S-196 

582.7 

2.000 

-93.4 

-.090 

-.430 

-.0108 

.3134 

18 

S-197 

565.9 

2.600 

-34.4 

-.150 

-.730 

-.0429 

.3543 

19 

S-158 

235.3 

-2.100 

8.0 

-.170 

.700 

-.0729 

-.0749 

20 

S-162 

581.2 

-1.700 

-29.1 

-.410 

-.070 

-.2190 

.4662 

21 

S-161 

190.5 

-5.800 

17.0 

-.460 

-1.080 

-.0774 

.0137 

22 

S-165 

376.8 

-2.300 

72.4 

-.720 

-1.360 

-.0345 

.1321 

23 

S-167 

147.2 

.400 

30.7 

-.570 

-.250 

-.0300 

.0690 

24 

S-172 

216.0 

-5.200 

-91.1 

-.130 

-.930 

-.0801 

.0344 

29 

S-169 

214.7 

-1.300 

-39.0 

-.130 

.310 

-.0192 

.0477 

26 

S-176 

588.3 

-10.100 

2.0 

-1.010 

. .020 

-.0233 

.1267 

27 

S-174 

169.9 

-3.100 

-38.2 

-.170 • 

-.690 

-.0696 

-.0159 

28 

S-175 

294.4 

-.700 

7.6 

-.070 

.110 

-.2135 

.0215 

29 

S-173 

115.5 

4.500 

-93.0 

-.080 

-1.050 

0.0000 

0.0000 

30 

S-177 

323.3 

.600 

-79.4 

-.630 

-.940 

.2290 

.0386 

31 

S-180 

486.0 

4.700 

-202.0 

-.240 

.060 

.2448 

-.1021 

32 

S-163 

120.0 

-.400 

-129.0 

.490 

-.120 

.1882 

.7061 

33 

S-183 

297.0 

-9.200 

47.0 

-.470 

0.000 

.0015 

.0170 

34 

S-184 

263.5 

6.700 

-8.0 

.260 

.320 

-.1355 

-.0692 

39 

S-182 

449.6 

8.000 

-177.0 

-.020 

.130 

.2266 

.1281 

36 

S-170 

299.7 

0.000 

-29.5 

-.760 

-.600 

.2933 

-.0426 

37 

S-189 

191.2 

-3.500 

70.5 

-1.140 

.340 

-.0743 

-.0923 

38 

S-181 

129.5 

-3.400 

84.2 

-.410 

-.280 

-.0777 

.0682 

39 

S-178 

604.4 

14.500 

-220.2 

-.030 

.170 

.3180 

.0590 

40 

S-190 

123.1 

2.400 

98.2 

.030 

.070 

-.0017 

.0943 

41 

S-188 

405.0 

-2.900 

149.5 

-.470 

-.620 

-.1845 

.1739 

42 

S-191 

219.2 

.600 

45.7 

-.960 

-.200 

-.0897 

-.2038 

43 

S-186 

124.2 

-4.200 

3.9 

-.240 

.690 

-.0147 

-.9603 

44 

S-189 

210.5 

-5.000 

61.1 

-.050 

.280 

.0101 

-.3837 

49 

S-187 

2 70.2 

4.000 

39.1 

.190 

0.000 

.0601 

.1696 

46 

S-194 

271.0 

.800 

14.1 

.030 

-.430 

.1489 

.0124 

47 

S-195 

194.8 

3.400 

-29.8 

-.090 

.770 

.0623 

-.1025 



SINPUTO 

TIUEI-TRAJECTORY to 200 NH ALTITUDE 
TITLE2<UITH PROPAGATION 900 SAMPLES 
Sl> 546679; 18989* 13. 9 >178.9* 30.69 «-120.a6« 
S2>1218031*17899* .5*179.0*16.48*>121.83» 

TC0&ST-300* 

I£RR0R-2*NORN-O« 

NSAHP-500* 

EhA6-l*3.6*3.6*0* 

NERPOP-3* 

SEHl- 20* 11 «6*>. 001 *-.001*. 009*. 010* 
SEN2«2550*>9.6* . .35* 0*. 001*. 002* 

SEN3- 1*«3.0* 0* .35*. 002*. 001* 

SEND 


DATE IS 06/30/81 
tine is 13.27.34 


AlO- 



routine step 

OTO NAR 1981 
PROS. NO. 2 

TRAJECTORY TO 200 NH ALTITUDE 
NO PROPAGATION 900 SAMPLES 


FLIGHTS FOR SCOUT FLIGHT EXPERIENCE ACCURACY 
STANDARD DEVIATIONS 

AIT-VEL-GAM-ZET-LAT-LONG <FT»FPS»DE6> 

29297. 80.4 .402 .925 .1641 .2239 

4.822 NN 

MEAN VALUES 

• ALT-VEL-GAM-ZET-LAT-LONC 

2391. -11.9 -.260 -.192 -.0347 .0427 


CORRELATION COEFFICIENTS 
ALT-VEL-6AH-ZET-LAT-L0NG 

-.43278 .43300 .11934 .40168 .14018 

-.32800 -.10826 -.33822 -.20899 

.16758 -.21837 -.11681 

.07191 -.47921 

-.07939 


COVARIANCE MATRIX lA MATRIX) 

ALT VEL GAM ZET LAT LONG 

.89830E-I-09 

-.10200E407 .64719E04 

.91043E404 -.10618E402 .16191E400 

.18399E404 -.49721E401 .39398E-01 .275588400 

.19316E404 -.44661E401 -.14423E-01 .61961E-02 .26943E-01 

.91966E403 -.37578E401 -.10529E-01 -.56335E-01 -.29183E-02 .90149E-01 


EIGENVALUES 

.85830E409 .52997E404 .28973E400 .13160E400 .35253E-01 .11113E-01 

EIGENVECTORS 

.10000E401 .11884E-02 -.21914E-05 -.39473E-05 -.11341E-06 -.31095E-05 

-.11884E-02 .10000E401 .44569E-03 .61192E-03 .39731E-03 .80723E-03 

.99470E-05 -.86939E-03 .16123E400 .95227E4|00 .40492E-02 .25918E400 

.21389E-09 -.49498E<^3 .95546E400 -.17894E400 .22700E400 .59516E-01 

.22905E-09 -.41268E-03 -.72590E-02 -.23903E400 -.39090E400 .88883E400 

.10719E-05 -.90665E-03 -.24707E400 -.63533E-01 .89199E400 .37319E400 

SENSITIVITY COEFFICIENTS 
ALT-VEL-GAM-ZET-LAT-LONC (TOP TO BOTTOM) 

.29297E409 .86188E-01 -.11799E-09 -.14319E-09 -.21294E-07 -.32780E-06 

-.34816E402 .72524E402 .23990E-03 .22198E-03 .74998E-04 .89096E-04 

.17423E400 -.62761E-01 .867-85E-01 .34545E400 .76026E-03 .27322E-01 

.62691E-01 -.32968E-0I .91429E400 -.64912E-01 .426Z2E-01 .62740E-02 

.65933E-01 -.Z9929E-01 -.39C73E-02 -.86710E-01 -.73394E-01 .93698E-01 

.31391E-01 -.36744E-01 -.13299E400 -.23047E-01 .16748E400 .39340E-01 
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SINPUTD 

TlTLEl-TRAJECTORY TO 200 NH ALTITUDE 
TITLE2-NQ PROPAGATION 900 SAMPLES 

TCOAST-0# • 

SEND 


DATE IS 06/30/81 
TINE IS 13.27.36 



ROUTINE STER 
OTO NAR 1981 
PROS. NO. 2 

TRAJECTORY TO 200 NM ALTITUDE 
NO PROPAGATION 500 SAMPLES 


FLIGHT RESULTS AT INPUT TINE 


SAMPLE 

VEHICLE 

PREOALT 

DELTA 

DELTA 

DELTA 

DELTA 

DELTA 

DELTA 

NUMBER 

NUMBER, 


ALT 

VEL 

GAM 

AZ 

LAT 

LONG 



N.NI. 

N.MI. 

FPS 

0E6 

DEG 

DEG 

DEG 

1 

S-136 

562.8 

9.700 

-58.5 

.770 

.280 

0.0000 

0.0000 

2 

S-131 

612.5 

6.900 

-280.0 

-.360 

-.150 

0.0000 

0.0000 

. 3 

S-138 

383.8 

3.200 

-107.9 

.150 

.770 

0.0000 

0.0000 

A 

S-139 

404.3 

8.900 

-29.2 

.050 

.190 

0.0000 

0.0000 

5 

S-140 

506.8 

-.200 

67.2 

-.370 

-.910 

0.0000 

0.0000 

6 

S-142 

487.1 

4.100 

142.9 

-.860 

-.290 

0.0000 

0.0000 

7 

S-143 

482.2 

4.200 

112.3 

.140 

-.280 

0.0000 

0.0000 

8 

S-145 

200.1 

-3.500 

15.0 

-.750 

.480 

0.0000 

0.0000 

9 

S-146 

495.2 

-3.700 

-24.3 

-.630 

0.000 

0.0000 

0.0000 

10 

S-147 

350.6 

3.800 

11.1 

-.980 

-.24a 

0.0000 

0.0000 

11 

S-148 

197.7 

2.800 

-31.4 

-.170 

-.550 

0.0000 

0.0000 

12 

S-149 

576.3 

0.000 

-38.6 

.160 

-1.160 

0.0000 

0.0000 

13 

S-150 

178.5 

1.000 

-31.1 

-.630 

-.140 

.0638 

.0339 

14 

S-154 

577.0 

-.700 

-57.5 

.120 

.260 

-.0473 

-.0361 

IS 

S-153 

117.1 

-.500 

-56.9 

-.630 

-.130 

0.0000 

0.0000 

16 

S-155 

279.4 

-5.50a 

94.3 

-.060 

.180 

-.1042 

.0115 

17 

S-156 

582.7 

2.000 

-53.4 

-.050 

-.430 

-.0108 

.3134 

18 

S-157 

565.9 

2.800 

-34.4 

-.150 

-.730 

-.0429 

.3543 

19 

S-158 

235.3 

-2.100 

8.0 

-.170 

.700 

-.0729 

-.0749 

20 

S-162 

581.2 

-1.700 

-29.1 

-.410 

-.070 

-.2190 

.4662 

21 

S-161 

190.5 

-5.800 

17.0 

-.460 

-1.080 

-.0774 

.0137 

22 

S-165 

376.8 

-2.300 

72.4 

-.720 

-1.360 

-.0345 

.1321 

23 

S-167 

147.2 

.400 

30.7 

-.570 

-.250 

-.0300 

.0690 

24 

S-172 

216.0 

-5.200 

-91.1 

-.130 

-.930 

-.0801 

.0344 

25 

S-169 

214.7 

-1.300 

-35.0 

-.130 

.310 

-.0152 

.0477 

26 

S-176 

588.3 

-10.100 

2.0 

-1.010 

.020 

-.0233 

.1267 

27 

S-17A 

169.9 

-3.100 

-38.2 

-.170 

-.690 

-.0656 

-.0159 

28 

S-175 

294.4 

-.700 

7.6 

-.070 

.110 

-.2135 

.0215 

29 

S-173 

115.5 

4.500 

-93.0 

-.080 

-1.050 

0.0000 

0.0000 

30 

S-177 

323.3 

.600 

-79.4 

-.630 

— . 540 

.2250 

.0386 

31 

S-180 

486.0 

4.700 

-202.0 

-.240 

.060 

• 2448 

-.1021 

32 

S-163 

120.0 

-.400 

-129.0 

.490 

-.120 

.1882 

.7061 

33 

S-183 

297.0 

-9.200 

47.0 

-.470 

0.000 

.0015 

.0170 

34 

S-184 

263.5 

6.700 

-8.0 

.260 

.320 

-.1355 

-.0692 

35 

S-182 

449.6 

8.000 

-177.0 

-.020 

.130 

.2266 

.1281 

36 

S-170 

299.7 

0.000 

-29.5 

-.760 

-.600 

.2533 

-.0426 

37 

S-185 

151.2 

-3.500 

70.5 

-1.140 

.340 

-.0743 

-.0523 

38 

S-181 

129.5 

-3.400 

64.2 

-.410 

-.280 

-.0777 

.0682 

39 

S-178 

604.4 

14.500 

-220.2 

-.030 

.170 

.3180 

.0550 

40 

S-190 

123.1 

2.400 

56.2 

.030 

.070 

-.0017 

.0943 

41 

S'188 

405.0 

-2.900 

149.5 

-.470 

-.620 

-.1845 

.1735 

A2 

S-191 

215.2 

.600 

45.7 

-.560 

-.200 

-.0857 

-.2038 

43 

S-186 

124.2 

-4.200 

3.9 

-.240 

.690 

-.0147 

-.5603 

44 . 

S-189 

210.5 

-5.000 

61.1 

-.050 

.280 

.0101 

-.3837 

45 

S-187 

270.2 

4.000 

39.1 

.150 

0.000 

.0601 

.1696 

46 

S-194 

271.0 

.600 

14.1 

.030 

-.430 

.1489 

.0124 

47 

S-195 

194.8 

3.400 

-25.8 

-.090 

.770 

.0623 

-.1025 
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APPENDIX B 

fortran code listings 


PROGRA»< STEP ( INPUT, 0UTPUT»TAPE1,TAPE8*T4PE3«INpuT»Tape6-CUTPUT) 
SCOUT TRAJECTORY ERROR PROPAGATION PROGRAM 

cov • covariance matrix of alt, vel, gam, zet,lat, long. 

INPUT IF IERRQR-4. 

DATAG ■ TABLE OF STANDARD DEVIATIONS OF PATH ANGLE FOR NORMALIZING. 

(.428,.428, .428,.430,.43A,.A42, .453 6N ) 

DATAH • TABLE OF STANDARD DEVIATIONS OF ALTITUDE FOR NORMALIZING. 

(3.35,4.85,5.80,6.43,6.90,7.20,7.30 BN) 

OATAV • ‘TABLE OF STANDARD DEVIATIONS OF VELOCITY FOR NORMALIZING. 

(75.6,64.6,68.8,91.3,92.5,93.0,93.0 BN) 

EMAG • VALUES FOR ERROR SOURCES USED TO CALCULATE COVARIANCE 
MATRIX OF A SPIN STABILIZED STAGE. ARRAY OF 4. INPUT 
WHEN PROPAGATION IS USED. VALUES ARE THE RATIO OF 
STANDARD DEVIATION DESIRED TO THE STANDARD DEVIATION 
USED TO CALCULATE THE SEN1-SEN4 DATA. 

HINJ • TABLE OF ALTITUDES FOR NORMALIZING. 

( IOC., 2C0., 300. ,400., 500., 600., 700. BN) 


lERROR • OPTION FOR INPUTTING DATA ERRORS IN FLIGHT DATA BASE. 
■1 INPUT flight RESULTS OF ALT,VEL,GAM. 

>2 INPUT FLIGHT RESULTS OF ALT, VEL, GAM, ZET, LAT, LONG. 
>3 INPUT STANDARD DEVIATIONS (SIG) AND CORRELATION 
COEFFICIENTS (RHO) OF ALT, VEL, GAM, ZET, LAT, LONG. 

>4 INPUT COVARIANCE MATRIX (COV) OF ALT, VEL, GAM, ZET, 
LAT, LONG. 


NERROR 

NORM 

NSAMP 

RHO 

SENl 

SEN2 

SEN3 

SENA 

SIG 

51 

52 

TCOAST 

TITLEl 

TITLE2 


NUMBER OF ERROR SOURCES OF THE SPIJ4 STABLIZED STAGE. 
MAXIMUM OF 4. (3 BN) 

NON-ZERO VALUE NORMALIZES ALT, VEL, GAM ERRORS OF THE 
flight data base to the altitude of si. (1 BN) 

NUMBER OF SAMPLES USED IN THE MONTE CARLO ANALYSES. 
(5000 BN) 

CORRELATION COEFFICIENTS IN ORDER OF ALT, VEL, GAM, ZET, 
LAT, LONG. ARRAY OF 15. INPUT WHEN IERROR-3. 
SENSITIVITY OF SPIN STABLIZED STAGE BURNOUT STATE 
PARAMETERS TO ONE SIGMA MOTOR PERFORMANCE ERROR SOURCE. 
UNITS ARE STATE PARAMETER UNITS PER SIGMA. ARRAY OF 6 
ORDER IS ALT, VEL, GAM, ZET, LAT, LONG. 

SAME AS SENl EXCEPT ERROR SOURCE IS PITCH TIPOFF. 

UNITS ARE STATE PARAMETER UNITS PER DEGREE. 

SAME AS SENl EXCEPT ERROR SOURCE IS YAW TIPOFF. 

UNITS ARE STATE PARAMETER UNITS PER DEGREE. 

SAME AS SENl EXCEPT ERROR SOURCE IS TIME OF SPIN 
STABLIZED STAGE IGNITION. UNITS ARE STATE PARAMETER 
UNITS PER SECOND. 

STANDARD DEVIATION OF ALT, VEL, GAM, ZET, LAT, LONG. 

INPUT WHEN IERPOR-3. 

NOMINAL state PARAMETERS AT STAGE BURNOUT. IF 
PROPAGATION IS SELECTED, SI IS STATE AT BURNOUT OF 
-NEXT TO LAST STAGE. ARRAY OF 6. ORDER IS ALT, VEL, 

GAM, ZET, LAT, LONG. 

NOMINAL STATE PARAMETERS AT LAST STAGE BURNOUT. 

INPUT IF TCOAiT IS NON-ZERO. ARRAY OF 6. ORDER IS 
ALT, VEL, GAM, ZET, LAT, LONG. 

nominal COAST TIME TO PROPAGATE COVARIANCE MATRIX. 

INPUT ZERO IF PROPAGATION IS NOT DESIRED. (0. BN) 

TITLE information PRINTED AT TOP OF EACH PAGE OF 
OUTPUT. 72 characters MAXIMUM. 
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c 


REAL 'me AKH, HE ANV# r^EANG# meaNZ *HEANL« HE ANL fi 
COHMON /BLK^/ yHU»SlGHA»6AUSS»GN(16)»XK,ZPAN* ICNT 
COMHON /COMIC/ STAlfc(6)» TCALL# KERR» OELT, T, TZERO 
COHHON /Olfi/ ZERO>OM6»TWO>RAO,Pl,TWOPl»6M,FT»RE»CQN»RQTATE 
OIHENSION SIG(6)/X(6)*0X(6)>RH0(19)* R(36)« RR(6*6)» W(6) 
OIHENSIOM lwnPC(100)» W0R0(11»100) 

OIHENSION TITLE1(8>, TITLE2(8J# BEST(6»4) 
dimension. SENK6), SEN2(6)« SEN3(6)« SEN4(6) 

DIMENSION MINJ(7)/0ATAH(7),0ATAG(7»»0ATAV(7I 
OIHENSION EPAR(6)» EXIC(6)* Sl(6)» S2(6)» SS(6)> EHA6I4) 
DIMENSION SUMI(6)r SUH2(6l* SUH3(6) 

DIMENSION SUHS0I(6)> SUHS02(6)* SUMS03(6) 

DIMENSION CPl(i5>. CP2(15), CP3U5) 

DIMENSION CMATl(2l)» CMAT2(211 

DIMENSION DHl(lOO)* DLKIOO)* DL6K100)* DVldOO)# OGKIOO) 
DIMENSION DZ1(1U0)> OH2(100)* OV2(100)t OG2(100)» 0Z2(1C0) 
DIMENSION D(21)f SC(6 a6)» C0V(21) 

EQUIVALENCE (BEST(1/1)*SENX(1>>* 

« (REST(l>2)>SEN2(l))f 

A (BEST(l*3)»SEN3(in# 

A (8EST(l»4)>SENA(l)) 

EQUIVALENCE (R ( 1 ) >RR ( 1* 1 ) ) 

EQUIVALENCE (SIGH »SIG( 1))*(SI6V *SI6( 2))»(SIG6 *SIG( 3)) 

A *(SIGZ tSIG( 4))t(SIGL *SIG( 5))«(SIGLG *SIG( 6)) 

A >(RHOHV >RHO( l))#(RHOHG »RHO( 2) ) * (.RHOHZ »RHO( 3)) 

A »(RHOHL »RHO( A } ) » { RHOHL G# RHO ( 5))»(RHOVG #RHO( 6)) 

A ,(RHOVZ »RHO( 7n»(RHQVL »RHP( 8)) , ( RH0VL6,RH0{ 9)) 

A *(RHOGZ *RH0(10) )» (RHOGL # RHO ( 11 ) ) « ( RHOGLG* RHO (12 ) > 

A »(RH07L fRH0(13))»(RH0ZLG>RH0(lA) )#(RHOLLG«RHO(15) ) 

DATA TITL£1»TITLE2/16*10H / 

DATA HINJ/UiC.>200«*30a.»AOb.f5CQ.>6CO.»700./* 

A 0ATAH/3t3 5»A.d‘>>5.e0#6.A57 6.90»7.2CW.30/> 

A 0ATAV/75.6»eA.A*88.8«91.3*92.5,93.0>93.0/> 

A OATAG/*A26».A28» .A28t .A30> .A3A> .AA2» .AS 3/ 

NAMELIST /INPUTO/ NSAMP ,NERROR» S I G, RHO» C OV .TCOAST. I ERROR , SI, S2» 
A SEN1,SEN2,SEN3,SENA,EMAG#HINJ,0ATAH,DATAG, 

A 0ATAV,N0RM 

C 

c 

REWIND 1 
REWIND 8 
CALL INIT 
NDATA-7 
NFLAG-3 
NSAHPaSOOO 
IERROR-2 
NORH-l 

TCOAST-0. 

MEAN-1 

NERROR-3 

NPRDB-0 

WRITE(6,330) 

C READ FIXED FIELD INPUT DATA 

CALL CAROS(N,100) 

RE AC (8, 3 20 > (IWORO(J), { WORD ( I, J ) , I -2, 11 ) » J *1# N ) 

10 CONTINUE 
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C READ NAHEIIST DATA 

CALL IMPUT(TITLfl»TlTLE2» 

REA0(8t INPUTO) 

NPR0B>NPRQ8-»1 

WRITE<6»3ACJ NPRQ8»TITLEX,TITLE2 
C INITIALIZE 

CALL INITl 
OQ 20 I<1*6 
SUH1(I)*0. 

SUN2(I)-0. 

SUM3(n-0. 

SUNS01(I)>0. 

SUNS02(I)>0. 

SUHSQ3(I)-0. 

20 CONTINUE 

DO 30 I-ltl5 
CPl(l)-0. 

CP2(I )>0. 

CP3(I)>0. 

30 CONTINUE 
C 

00 AO I-l»6 
STATEiD-SKI) 

AO CONTINUE 
CALL CONIC 

0ELT-ANAXl(TC0AST»l.E-3) 

CALL TSTEP 
00 SO I-l>6 
SS(I)«STATE(I) 

50 CONTINUE 
C 

ALT-S1(1)/FT 

ZET-SKA) 

XLAT-SKS) 

C 

C SET-UP OP INPUT DATA ACCORDING TO lERROR 

GOTO (60»60>160«170)> lERRCR 
C 

60 CONTINUE 

SZET-SIN(ZET/RA0) 

CZET«COS(ZET/RAO) 

CLAT-COS(XLAT/RAO) 

WRITE (6« 350) 

WRITE (6* 360) ( J * IWORO ( J ) » ( WORD ( I » J ) > I >3» 9 ) # J >1» N ) 
WRITE(6,3A0) NPRQ8»TITLEi,TlTLE2 
WRITE(6»370) ALT# ZET.tLAT 
WRITE(6#380) 

CALL INTER ( ALT#0EVH1#NDATA#0ATAH»HINJ ) 

CALL INTER ( ALr>OEVVl>NDATA#DATAV#HINJ ) 

CALL INTER (ALT#0EVG1#N0ATA*0ATAG*HINJ ) 

Nl-0 

N2-0 

N3«0 

C 

DO 120 I-lfN 
XH-WORO( A,I)*FT 
XV>WORO( 5# I ) 
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X&-W0R0'(6* I ) 

XZ-W080(7.I) 

XL>VOPO(8»n 
XLG«wn80(9* 1) 

IF (NQRH.EO.O) GHTH 70 . 

c nqrhalize flight errors 

XALT«W0R0(3»I) 

CALL INTER ( XALT»0EVH2»N0ATA»0ATAH,HINJ ) 

CALL INTER (XALT»0EVV2»N0ATA»0ATAV»HINJ ) 

CALL INTER ( X ALT « 0EVG2* NOAT A» 0ATA6> HIN J ) 

XH«0EVHl/0EVH2*XH 
XV»0EVV1/0£VV2*XV 
XG-0EVGI/DEVG2AXG 
70 IC00E-WQR0(2>I) 

GOTO (80«80«90»I20)* ICOOE 
C 

C CONVERT LAT^LONG DEVIATIONS TO RANGE rCROSSRANGE DEVIATIONS 

80 SZ>SIN(UOROai#l) /RAO) 

CZ>C0S(W0R0(11,I)/RA0) 

cn>re«xl/rao 

CE>RE*XLG*C0S(W0R0<10»I)/RAD)/RAD 
OR»(CN*CZ+CE*SZ) /FT 
OCR-(CN*SZ-CE*CZ)/FT 

CONVERT RANGE#CROSSRANGE DEV TO LAT>LONG DEV FOR INPUT 
AZINUTH AND LATITUDE 
CN-0R*CZET+DCR*5ZET 
CE-OP*SZET-DCP*CZET 
XL»CN/RE*FT*RAD 
XL6-CE/RE*FT/CLAT*RA0 
90 CONTINUE 
N3-N3+1 

WRITE(6*390) N3» XHf XV» XG# XZ » XL> XLG> ICOOE 
GOTH (1C0#100»110*120) * ICOOE 

lOU Nl-Nl+l 

PHKNl) «XH 
0V1(N1)-XV 
DG1(N1)«XG 
DZKND-XZ 
0L1(N1)-XL 
OlGKNl) -XLG 

GOTO a20*110»110«12U)» ICOOE 
110 N2-N2+1 

0H2(N2) -XH 
0V2<N2)«XV 
0G2(N2)«XG 
0Z2{N2)«XZ 
120 CONTINUE 

WRITE(6,3<.0) NPR0BtTITLEl»TITLE2 
IF (lERRQR.EO.l) GOTO 1^0 
IF (lERROR.GT.l) GOTO 150* 

130 CALLCORCO ( DH2» 0 Z2 » N2» S IGH» SIGZ» RHOHZ» HE ANH» NEANZ» HE AN ) 

CALL CORCn (OV 2 »OZ 2 »N 2 »SIGV»SI 6 Z/RHOVZ»HEANV»HEANZ#HEAN) 

CALL CORCO (D62»0Z2»N2,SIGG»SIGZ#RH0GZ»>*EANG»HEANZ»HEAN) 

NFlAG-0 

WRITE(6f3A0) NPRQBtTITLEl#TITLE2 
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140 WR1TE(6»400) 

CALL CORCn (OH2,OV2>N2>S1GH*SIGV>RHOHV»HEANH,HEANV*HEAN> 

CALL CORCO (0H2,0G2»N2,SlGH»SI66>RH0HG»HEANH«>1EANG>r«EAN) 

CALL CQRCQ ( 0V2>0G2 *N2, SIGV» SIGG» RHOVG* HEANV » HE ANG> f«E AN ) 

GOTO 160 

ISO CALL CQRCQ < 0H1> OVlf Nl# S 1GH> S IGVyRHPHV# HE ANH* HE ANV* HE AN ) 

CALL CORCQ ( OHl • 0G1*N1« SIGH» SIGG* RHOHG# HE ANH, HE ANG* HE AN ) 

CALL CORCQ ( nHl,0 Z1«NI> SIGHTS I6Z» RHQHZ* HE ANH » HE ANZ » MEAN ) 

CALL CORCQ ( 0H1»0L 1 *N1 » SIGH>S IGL r RHOHL# HE AKH, ME ANL » HE AN ) 

CALL COPCO (OHl*DLGl«Nlf SI6H»SIGL6>PH0HLG»HEANH»HEANLG»HEAN) 
CALL CORCO (OV1,OG1»N1»SIGV»SIG6 *RHOVG>heanV»HEANG»HEAN) 

CALL CORCO <DVl»nzi«Nl,SIGV>SIGZ^RMOVZ«HEANV»HEANZ»HEAN) 

CALL CORCO (0V1»0L1>N1*SI6V#SIGL>RH0VL»HEANV»HEANL>KEAN) 

CALL CORCO (OVlfDLGl»Nl,SIGV>SIGLG»RHOVL6»HEANV»HEANLGtHEAN) 
CALL CORCO (OG1>OZ1«N1>SIG6#SIGZ>RHOGZ»HEANG»HEANZ»HEAN) 

CALL CORCO (OGl»OLi*Ni«SIGG>SIGL»RHOGL»HEANG>HEANL»HE^AN) 

CALL CORCO (CGlf 0LGl#Nl,SIG6f SIGLG»RHOGLG»HEANGfHEANLG>HEAN) 
CALL CORCO <071*OL1*N1*SIGZ#SIGL»RHOZL>HEANZ«HEANL*HEAN) 

CALL CORCQ ( CZl> OLGlt N1 » SIGZ» SIGLG* RH0ZLG»HE ANZ» HE ANLG* HE AN ) 
CALL CORCO (OLl>OLGl»Nl*SIGLfSlGLG>RHOLLG«HEANL*HEANLGf HEAN) 
NFLAG-1 
WRITE(6#410) 

WRITE (6*420) N1 
160 D(1)-SIGH«SIGH 

D(2)-RH0HV*SIGH*SIGV 

0(3)»SIGV4SI6V 

0(4)«RH0HG*SIGH*SIGG 

0(5)»RH0VG*SIGV*SIGG 

0(6)>SIGG*SIGG 

0(7)»RH0HZ*SIGH4SIGZ 

0(8)«RH0VZ*SIGV*SIGZ 

0(9)«RH0GZ*SIGG*SIGZ 

D(IO) -SIGZ4SIGZ 

0(11)«RH0HL*SIGH*SIGL 

0(12) •RHOVL*SIGV*SIGL 

D(13)*RH0GL*SIGG*SIGL 

O(14)-RHOZL*SIGZ*SIGL 

D(15)-SIGL4S1GL 

0(16) -RHOHLG*SIGH*SIGLG 

0(17)«RHQVLG*SIGV*SIGLG 

0(18)-RH0GLG«SIGG«S1GLG 

0(19)-RN0ZL64S1GZ*SIGLG 

0(2C)-RH0LLG«SIGL4SIGLG 

0(21)>SIGLG*S16LG 

WRITE (6*430) S I GH*S IGV* S IGG*S IGZ*S IGL* S IGLG 
SI GALT- SIGH/60 76. 11549 
WPITE(6*440) SIGALT 

WRITE (6* 450) HEANH* HEANV*HEANG»HE ANZ * HEANL *HE ANLG 
WPITE(6*460) RHOHV* PHOHG* RHOHZ* RHOHL* RHOHLG* 

* RHOVG* RHOVZ* RHOVL* RHOVLG* 

A RHOGZ* RHOGL* RHOGLG* 

A RHOZL* RHOZLG* 

A RHOLLG 

GOTO 190 

170 00 IPO 1-1*21 
leO 0(I)-C0V(I) 

190 WPITE(6*470) 
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WRITE(6»4B0) D 
IF (NFLAG.FO.i) GOTO 13C 
CALL EIGEN <0>R»6>0) 

WRITE(6» AQO) 

C 0 ARRAY IS NOW THE EIGENVALUES 

C RR ARRAY IS THE EIGENVECTORS 

WRITE (6* 500) D(1)»00)>D(6)*0(10)*0(15)»0(21) 

WRJTE(6>510) 

WRITE (6* $00) ((RR(I*J)»J-1»6)»I*1«6) 

Wm«SQRT{0(l)) 

W(2)«S0RT(0(3) ) 

V( 3)"S0RT(0(6) > 

W( A)«SQRT(D (10) ) 

V(5)«SQRT(0(1S) ) ' 

■ W(6)«SQRT(D(2I) ) 

00 200 

no 20C K>1,6 

SC (J#R)«W(K)*RR(J,K) 

200 CONTINUE 

WRITf(6,520) 

WRITE(6»500) ( (SC(I>J)»J-1>6)#I«1*6) 

C 

IF (TCOAST.EO.O.) GOTO 310 
C 

WRITE(‘>»3A0) NPR0B»TITLE1»TITLE2 

C • COMPUTE LAST STAGE BOOST COVARIANCE MATRIX CMATl 

00 230 NTIMES-1»NSAMP 
CALL NQRRAN 
X(1)«GAUSS 
CALL NORPAN 
T0MAG>A6S(GAUSS) 

CALL UNIPAN 

thet«zran*twopi 

X(2)-TOMAG*COS(THET) 

X(3)«T0MAG*SIN(THET) 

IF (NERRCR.LT.A) GOTO 210 
CALL NORRAN 
X(A)>GAUSS 
210 CONTINUE 

C DX IS LAST STAGE BOOST STATE ERRORS 

00 220 I>1»6 
DX(I)>0. 

00 22C J«1»NERR0R 

220 OX (I )«SEST(1»J )*X ( J )*EMA6( J)+OX( I) 

CALL COVR (CMATl, OX, NTIMES»SUM1»SUMS01#CP1) 

230 CONTINUE 
C 

CALL COVRl (CMATl, OX, NSAMP,SUM1,SUMS01, CPI) 

WRITE(6,560) 

WRITF(6,$30) 

WRITE(6,AB0) CMATl 


COMPUTE COVARIANCE MATRIX AFTER PROPAGATION 
00 2P0 NTIMES-1»NSAMP 
C sample input epoch ERRORS 
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00 2<»0 I » l *6 
CALL NORRAN 
2A0 EPAR(I)-6AUSS 
00 2S0 J>1#6 
FXIC( J)-ZER0 
00 250 I>1*6 

250 EXIC(J)-EXIC(J)«SC(J>n*ERAR(I) 

00 260 I>1*6 

C STATE IS AT INPUT EPOCH 

260 STATE(I)-Sl(n*EXIC(I) 

CALL COVR (0»EXIC,NTIMES*SUH2,SUMS02»CP2) 

C INITALIZS ORBIT PARAMETERS FOR STATE 

CALL CONIC 
OELT-TCOAST 

•C ADVANCE STATE 

CALL TSTEP 

C OX IS state ERRORS AFTER PROPAGATION 

DO 270 I-l»6 

270 OX(I)-STATE(I)-SS(n 

CALL COVR (CMAT2#DX,NTIMES»SUM3»SUHS03»CP3) 

280 CONTINUE 
C 

WRITE(6f 5A0) 

WRITE(6*530) 

CALL COVRl (0 *0X»NSAMP»SUM2>SUMSQ2*CP2) 

VRITE(6«A80) 0 

WRITE(6>550) 

WRITE(6>530) 

CALL COVRl <CMAT2>0X»NSAMP,SUM3*SUMS03»CP3) 
WPITE(6fA80) CMAT2 


COMPUTE COVARIANCE MATRIX AT LAST STAGE BURNOUT 

00 290 I-l#21 
290 D«I)-CMAT1(I)+CMAT2(I) 

WRITE(6»3A0) NPROB# TITLEl# TITLE2 
WR1TE(6>570) 

WRITE(6f 530) 

WRITE(6>480) 0 


COMPUTE SENSITIVITIES 


CALL 

EIGEN ( 

0* 

WR 

ITE 

(6*490) 


UR 

ITE 

(6*500) 

0 

UP 

ITE 

(6* 510) 


UR 

ITE 

(6*500) 

( 

U( 

1)- 

S0RT(0( 

1) 

U( 

2)- 

SORTLO( 

3) 

U( 

3)- 

S3RT(0( 

6) 

U( 

A)« 

SQRT(0( 

10 

W( 

5)- 

SORT(0( 

15 

U( 

6>» 

SORT(0( 

21 

00 

300 J-l*6 


DO 

300 K«l*6 



10)#0(15)/0(21 ) 
I-l»6) 
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300 SC(J»K)>W(K)*RR(J«K) 

WQITE(6*520) 

WRITE (6« 500) ( (SC (I*J )«J-1*6)*I>1»6) 

C 

310 CONTINUE 

C WRITE SENSITIVITIES ON UNIT 1 FOR READING BY SOAR 

WRITE (1) ( (SC(I*J)«J-1»6)»I-1«6) 

GOTO 10 


320 FORNAT ( IX, AA,FA.0» 5F7.0»2F8.0*2F10.0) 

330 FORMAT (♦!♦) 

3A0 FORMAT (♦!♦, T72, ♦ROUTINE STEP*/T72» ♦DTD MAR 1981P/ 

* T72,*PR08. NO. *13/ 10X,8Ai0/ 10X,8A1C//) 

350 FORMAT (28X,*FLIGHT RESULTS AT INPUT TIME*//6X, ♦SAMPLE VEHICLE P 
• aRE0ALT*,5X, ♦DELTA DELTA DELTA DELTA DELTA DELTA*/6X, *NUH 

aBER number*, 15X,*ALT VEL gam A2 LAT LONG*/, 2 

aAX,*N,MI.*,6X,*N.MI. FPS*,5X,*0EG*,5X,*0E6*,5X,*DEG*,5X,*DE6*/) 

360 FORMAT ( 6X, 13, 6X , *S*, A^, F9 .1, Fll .3, F8 .1, 2F8 . 3, 2F8 .4 ) 

370 FORMAT ( 12X*N0RMAL I ZED DEVIATIONS AND ADJUSTED LATITUDE, LONGITUDE 
A OEVIATIONS*/20X,*ALTITUOE"*,T38,F9.1,* N MI*/28X,* AZ IMUTH-*, T38, F 
a9.3,* DEG*/28X,*LATITUDE«*,T38,F9.4,* DEG*//) 

380 FOR**AT (7X, ♦SAMPLE ALT VEL GAMMA ZETA LAT 

A LONG C0DE*/7X, ♦NUMBER FT FPS DEG D 

aEG DEG DEG*/) 

390 FORMAT! Ill, F12.0,F10.1,2F10.3»2F9.4, 18) 

400 format (20X, ♦flights FOR SCOUT FLIGHT EXPERIENCE ACCURACY*/) 

410 FORMAT ( 20X, *FL IGMTS WITH LATITIUDE /LONGITUDE DEVIATIONS*) 

420 format (30X,12»* FLIGHT SAMPLES*) 

430 FORMAT ( lOX, *STANOARO 0EVIATIPNS*/10X*ALT-VEL-GAM-ZET-LAT-L0NG (FT 
a,FPS»OEG )*//10X, F10.0,F10.1,2F10.3,2F10.4) 

440 FORMAT (10X,F10.3,* NH*/) 

450 format (10X/*MEAN VALUES*/10X,*ALT-VEL-GAM-ZET-LAT-L0NG*/ 

A 10X,F1C.0,F10.1,2F10.3,2F10.4//) 

460 FORMAT ( lOX^CORRELATIQN COEFF IC lENTS* / lOX* ALT-VEL-GAM-Z E T-L AT-L0N6 
a*//10X,5F10.5/20X,4F10.5/30X,3F10.5/4uX,2F10.5/50X,F10.5/) 

470 format (/25X*C0VARIANCE MATRIX (A MATR IX ) */12X*ALT*, 9X*VEL*, 9X, *GA 
aH*,9X*ZET*,9X*LAT*,9X*L0N6*/I 

4 80 FORMT ( 6X,E13.5/6X,2E13.5/6X,3E13. 5/6X,4E13. 5/6X,5E13.5/,6X,6E13. 
a5// ) 

490 FORMAT ( 31X , *E IGENV ALUES* ) 

500 format (6X,6F13.5) 

510 FORMAT ( /31X,*EIGfcNVECTQfiS*/) 

520 FORMAT ( /25X*SENSITIVITY CaEFFICIENTS*/25X*ALT-VEL-6AM-ZET-LAT-L0N 
aG (TOO TO BOTTOM)*/) 

530 format (13X*ALT*10X*VEL*10X*GA'1*10X*Z£T*10X*LAT*9X4LONG*) 

540 format (//1uX,*SAMPLE0 COVARIANCE MATRIX AT INPUT EPOCH*/) 

550 FORMAT ( //10X,*C0VARIANCE MATRIX AFTER PROPAGATION*/) 

560 FORMAT ( //10X,*C0VARIANCE MATRIX OF LAST STAGE BOOST*/) 

570 format (//10X,*COVARIANCE MATRIX AT LAST STAGE BURNOUT*/) 

END 


- B8 - 



o o o o 


subroutine CAP0S(NSANP>NHA«) 

THIS subroutine READS FIXED FIELD DATA CARDS IN (A) 

F0RNAT» WRITES ON UNIT 6 FOR PRINTOUT AND WRITES ON 
UNIT 8 FOR SUBSEQUENT READING IN (F) FORNAT. 

IMPLICIT INTE6EH(A-Z) 

DIMENSION CARDtB) 

NSAMP-D 

MAX-NNAX*! 

REA0(5#50) CARO 

IF (EOF(S) .NE.O) GOTO 30 
WRITE(6»60) CARO 
DO ID I>1»MAX 

REA0(5«S0) CARO 
WRITE(6*60) CARO 
WRITE(6«50) CARD 
IF (CARO(l) .EQ.AHENO ) GOTO 20 
10 CONTINUE 

URITE(6*70) NNAX 
GOTO AO 
20 CONTINUE 
NSAMP-I-1 

IF (NSAHP.LE.O) WRITE(6»801 
IF (NSAMP.LE.O) GOTO 40 
REWIND 8 ^ 

RETURN 
30 CONTINUE 
WRITE(6«90) 

40 STOP 

90 FORMAT (8A10) 

60 FORMAT (5X*PA10) 

70 FORMAT ( //15X, *SU8R0UT IN6 CARDS - SAMPLES EXCEEDS MAX OF ♦ 14) 
80 FORMAT ( //l9X,*SL'BRnUTINE CAROS - NO DATA PROVIDED*) 

90 FORMAT < //15X»*SUBR0UTINE CAROS - NO SAMPLES INPUT*) 

END 
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SUBBOUTtNE CONIC 

INITIALIZES A CONIC FROM A SPHERICAL STATE VECTOR (STATE) 
FOR SUBSEQUENT PROPAGATION OF THE STATE VECTOR ALONG THE 
CONIC* WHICH IS ACCOMPLISHED BY SUR. TSTEP. 

COMMON /CENT/ C6 TA* SETA ,CG* S6»R *V* SEA 
common /CONI.C/ STATE(6»* TCALL* KERR* STEP* T* TZERO 
COMMON /CORO/ 0(3*3) 

COMMON /DIG/ ZER0*QNE,TW0*RA0*PI*TU0PI*GM»FT»RE*C0N, ROTATE 
common /man/ A,A1EE* AVG*E*006*6AREF»ETA» icircl*lousy*p* 
PGME*RREF»TAN0N*TRUE»TREF*XNA*PERI00,6P 
DATA SMALL/1. E-6/ 

c 

T-TZERO 
ICIRCL - 0 
LOUSY > 0 

C 

SG • SIN( STATEO) * CON ) 

CG - COS( STATEO) ♦ CON ) 

SZ • SIN( STATE(<») ♦ CON ) 

CZ • COS( STATF(A) ♦ CON ) 

SL • SIN( STATE(5) * COM > 

CL • COS( STATE(5) ♦ CON ) 

TEMP • STATE(6) ♦ CON 
SO • SIN( TEMP ) 

CO • CnS( TEMP ) 

P • STATEd) +RE 
V ■ STATc(2) 

c 

0(1*1) > CO * CL 

0(1,2) • - SO ♦ SZ - CO • SL ♦ CZ 

0(1*3) • SO ♦ CZ - CO • SL * SZ 

0(2*1) • SO ♦ CL 

0(2*2) • CO * SZ - SO ♦ SL ♦ CZ 

0(2,3) - - CO * C2 - SO ♦ SL • SZ 

0(3*1) • SL 

0(3*2) ■ CL * CZ 

0(3*3) • CL « SZ 

C 

TP • ZERO 
TA • ZERO 
C 

A > R * 6M / ( TWO ♦GH-R*V*V) 

IF (A.GT.C.) GOTO 30 
10 LOUSY • 1 

C COMPUTED A NON-ELLIPTIC* OR LOUSY* CONIC. 

C IT WILL BE IGNORED. 

WRITE (6*20) 

20 FORMAT ( /* 6X 20HL0USY CONIC FOR BODY ) 

RETURN 

30 OOG ■ R ♦ V • CG 
P - OOG ♦* 2 / GM 

E • SQRT( ONE - P / A ) 

IF (E-SMALL. GT. C. ) GOTO AO 
E > ZERO 

C CIRCULAR orbit ENCOUNTERED. SET ICIRCL TO UNITY. 
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ICIRCL - 1 



GOTO 50 

40 IF (ABS( F-ONE)-S*<ALL .LE. C. ) GOTO 10 
PG«E • SORT( P / GM } / e 
50 AVG • SORT! GH / A**3 ) 

PERIOD ■ TWOPI / AVG 
XNA ■ AVG ♦ A 

AIEE • A ♦ SORT( ONE - E**2 I 
IF UCIRCL.GT.O) GOTO 6C 
STA • PGME • V • S6 
CTA »(P-R )/(R *E) 

TA • ATAN2C STA, CTA ) 

SEA • R * STA / AIEE 
CEA • E ♦ R * CTA / A 
EA - ATAN2( SEA, CEA ) 

IF (EA.LT.ZERO) EA-EA>TUQPI 
TP • C EA - E * SEA J / AVG 
60 TANON • TA 
TRUE • TA 
RREF > R 
EAREF • EA 
TREF • TP 
ETA - ZERO 
RETURN 
END 


- Bll - 



SUni»nuTfN6 COPCG ( X , Y ,NS AMP ,$ IGX, S IGt » BHOlf Y , X MEAN, YME AN» MEAN ) 
CALCULATES STANOARO DEVIATIONS, mean VALUES AND THE 
CORRELATION COEFFICIENT OF TWO VARIABLES OF (NSAMP) 
RANDOM samples EACH. 

dimension xa),Y(l). 


XSUM-C. 

YSUM-0. 

xsosuM-a. 

YSOSUM-0. 

XYSUM-0. 

SAMPN-NSAHP 

DO lu I-l»NSAMP , 

XSUM-XSUM+Xd) 

YSUM-YSUM*Y(H 

xscsuM-xsosur+x { i )*x (i ) 

YSOSUM«YSOSUM^Y(I)*Y(I) 

XYSUM-XYSUM4X( I»*Y( I) 

XMEAN-XSUM/SAMPN 
YMEAN.YSUM/SAMPN 
IF (MEAN.EO.C) XMEAN-0. 

IF (MEAN.EQ.O) YMEAN*0. 

SIGX-SOBT (XSOSUM/ (SAMPN-l.)-XMEANPXMEAN) 

SI GY- SORT (YSOSUM/(SAMPN-l,)-YMEAN*YMEAN) 
RHnXY-(XYSUM-SAMPN*XMEAN*YMEAN)/SAHPN/SlGX/S IGY 
RETURN 
ENO 
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SUBPOUTINE CQVR(COV,OX,NTIHES»SUM*»SUMSay ,CP) 

GENERATES A SY»*METRIC COVARIANCE NATRI* OF SIX PARAMETERS. 
COVR IS CALLED EACH SAMPLE. 

COVRl IS CALLED AFTER ALL SAMPLES APE CALCULATED. 
COVARIANCE MATRIX IS GENERATED. 

DIMENSION 0X(6)*SUHX(6)>SUMS0X(6).CP(15)»XM0M2(6)» 

A RH0(15)tSI6(6).SVAR(6)*CQV(21) 

CALL WHEN EACH SAMPLE SET IS GENERATED. 

XN • FLOATINTIMESJ 
XNP • (XN - l.)/XN 
DO 10 I-1.6 
02X - OX(I)*OX(n 
SUMXdJ • XNP*SUMX(I) ♦ 0X(1)/XN 
SUMSOX(I) • XNP*SUMSOX(I J ♦ 02X/XN 
10 CONTINUE 
DO 20*I>lf9 
LB • !♦! 

DO 20 J>LB«6 
PROD ■ OX(I)*DX(J) 

INO • I ♦ (J*J -3*J ♦2)/2 
CP(INO) • XNP*CP(IND) ♦ PROD/XN 
20 CONTINUE 
RETURN 


CALL WHEN ALL SAMPLE SETS ARE GEN-ERATED AND 
COVARIANCE MATRIX IS DESIRED. 

ENTRY COVRl 
XN-FLOAT(NTIMES) 

XNPR • XN/(XN-1.) 

DO 30 I-1.6 

XM2 < SUMX(U*SUMX(I) 

XM0M2CI) - SUMSOX(I) - XM2 
SVAP(I) ■ XNPR*XM0M2(I) 

30 CONTINUE 
DO AO I>1»9 
LB • !♦! 

DO AO J-LB.6 

INO - I ♦ (J*J-3*J*2»/2 

DENOH - SORT(XMQH2(1)«XHON2(J) ) 

RHO(INO) • (CP(INO)-SUMX(I)*SUMX(J) )/OENQM 



IF (OENOM.EO.O. ) RHO 

(INO)-O. 

AO 

CONTINUE 




DO 90 I> 

1#6 


50 

SIG(I) > 

SORT(SVAR( 

in 


CQV(l) - 

SVAR(l) 



CQV(2) - 

RH0(1)«SI6( 

1)*SIG(2I 


C0V(3) > 

SVAR (2) 



COV(A) • 

RH0(2)«SIG( 

1J*S1G(3I 


C0V(5) • 

RH0(3)*SIG( 

2)*SI6(3) 


C0V(6) » 

SVAR(3) 



C0V(7) - 

RHO(A)«SIG( 

1)*SIG(A) 


C0V(8) • 

RH0(9»*SIG( 

2)*SIG(A) 


C0V(9)« 

RH0(6)*SIG(3 

)*SIG(A) 


COV(IO) 

• SVAP(A) 



COV(ll) 

• RH0(7)*SIG 

(1)*SIG(9) 
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CnV(12)*- RH0t8)*SIG(2)*SIG(5) 
C0V(13) - RH0(9)«SIG(3MSIG( 
C0V(14> • RH0(1C)*SIG(4J«SIG(5) 
C0V(1») • SVAR(5) 

COV(lft) • RH0(11)*SIG(I)*SIG(6) 
CQV(17) ■ RHO(l2)*SI6(2 )*SIG(6) 

rOV(l8) • RHO(13)*SIG(3)«SIG(6) 

C0V(19) - RH0(1«)*SIG(4)«SIG(6) 

C0V(20) • RH0(15>*SIG(5)*SIG<6> 

C0V(21) • SVAR(61 

RETURN 

ENO 
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SUBROUTINE EIGEN(i,R«N«HV) 
DIMENSION A(1)>R (1) 


DESCRIPTION OF PARAMETERS 

A - original symmetric MATRIX, DESTROYED IN COMPUTATION. 

RESULTANT EIGENVALUES ARE DEVELOPED IN DIAGONAL OF MATRIX 

A IN DESCENDING ORDER. 

R - RESULTANT MATRIX OF EIGENVECTORS STORED COLUMNWISE IN 

SAME .SEQUENCE AS EIGENVALUES. 

N * ORDER OF MATRICES A AND R 

MV - INPUT CODE 

0 COMPUTE EIGENVALUES AND EIGENVECTORS 

1 COMPUTE ONLY EIGENVALUES 

METHOD - DIAGONALIZATION METHOD ORIGINATED BY JACOBI AND ADAPTED 
BY VON NEUMANN FOR LARGE COMPUTERS AS FOUND IN MATHEMATICAL 
METHODS FOR DIGITAL COMPUTERS, EDITED BY A. RALSTON AND H. S. 
MILF, JOHN WILEY AND SONS, NEW YORK, 19»i2, CHAPTER 7. 

GENERATE IDENTITY MATRIX 

5 RANGE-1. OE-12 

IF (HV-1) 10,25,10 
10 10— N 

DO 20 J-1,N 
I0-I04N 
DO 20 I-1,N 
IJ-IQ>I 
RIIJ) - 0.0 
IF (I-J) 20,15,20 
15 R(lJ) • 1.0 
20 CONTINUE 

COMPUTE INITIAL AND FINAL NORMS ( ANORM AND ANORMX) 

25 ANORM-0.0 
DO 35 1-1,N 
DO 35 J-I,N 
IF <I-J) 30,35,30 
30 IA-I+(J*J-J y/2 

ANORM-ANORM-fAdA )«A(IA) 

35 CONTINUE 

IF (ANORM) 165,165,40 
40 ANQRM-1.414PSQRT (ANORM) 

ANRMX - ANORM*RANGE/ FLOAT(N) 

INITIALIZE INDICATORS AND COMPUTE THRESHOLD, THR 

IND-0 

THR-ANQRM 

45 THR • THR/FLOAT(N) 

50 1-1 
55 M-L*l 

COMPUTE SI.N AND COS 

60 M0-(M*M-M)/2 
L0-(L*L-L)/2 
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C 

C 

c 


lH«L+*'0 

<S2 IP (ABS (A(LM) )-THR )130»65»65 
65 INO-I 
LL-L+LO 
NM-n^MQ 

X-C.5*{A(LU-A(HHn 
6B Y--A(LM)/SORT ( A ( L»n ♦ A ( LM ) +X** ) 

IF tX) 70»75»75 

75 SINX«Y/SQ8T ( 2.0* ( l.O^ ( SOPT (1»0-Y*Y)))I 
S1NX2«5INX*SINX 

78 CQSX«$aRT (1,0-SINX2) 

CQSX2«CaSX^COSX 

SINCS«SINX*COSX , 

ROTATE L AND N COLUMNS 


1LQ>N*(L-1) 

IHO«N*{M-l) 

00 125 I"1#N 

IF (I-U 80»115,8« 

80 IF (I-M) 85>115»90 

89 IH-I+MO 
GOTO 05 

90 JN-H+ia 

95 IF (I-L) 100»105»105 
100 IL-I+LO 
GOTO llu 


lv.5 IL-L + IQ 

110 X-A(IU*C0SX-A{1M)*SIMX 

A(IH)-A( IL)*SINX*A{If*)^CDSX 


A(IU«X 

115 IF (MV-1) 120fl25»120 
120 ILR-UO+I 
IMR-IMO+I 

X-R(ILR)*COSX-R(IMR>*SINX 
R( IHP)-B( ILP)*SINX*R(IMR»*C0SX 
R(ILR )»X 


125 CONTINUE 

X»2.0*A{LM)*SINCS 

Y-A(LL)*C0SX2+A(MM)*SINX2-X 


X»A( LU *SINX2 + A ( MM) *C0SX2^X 

AILM) ■( A(LL »-A(MM) )*SINCS + AILH)*<C0SX2-SINX2) 
A(LL)-Y 


AO-Mj-X 


TESTS FOR COMPLETION 
TEST FOR N • LAST COLUMN 

130 IF (M-N) 135»1A0»135 
135 M-M+1 
GOTO 60 

lAO IF <L-(N-D) 1A5»150»145 
165 L-L+l 
GOTO 55 

150 IF (INO-1) 160*I55#160 
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155 INO-0 
GOTO 50 


COHPARfe THRESHOLD WITH FINAL NQRM 

160 IF ITHR-ANRHX) 165»L6S#A5 
C SORT EIGEN VALUES AND EIGENVECTORS 

165 10— N 

00 185 I*1>N 
lO-IO+N • 

u-uci*i-n/2 

J0«N*(I-2) 

00 185 J«I»N 

FN«J*U*J-J)/2 

IF fA(LL)>A(HH)) 170#185»185 
170 X-ACLU 

A(LL)-A(NH) 

A(HM).x 

IF (NV-1) 175»185»175 
175 DO 180 K>1>N 
ILR>I0-^K 
IHR-JQ+K 
X"R(ILR) 

R(ILR)-R(INR) 

180 R(IHR)-X 
185 CONTINUE 
RETURN 
END 
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SUBROUTINE INIT 
INITIALIZES CONSTANTS 

COHNON /DIG/ ZERa«ONE,Twa>ftAO»PI>TUOPI>GN,ET^RE«CON»ROTATE 
COHHQN /BLK4/ XHU* SI GNA , GAUS S>6N( 16 ) > X K , ZR AN , ICNT 
DATA GN/23.*71.»47.»A3.»S9.>19.*61.*37.>89.«97.#13.#29.*73.«83* 
A 31. #67. / 


ZERO-0.0 
ONE*l .0 
Twn-2.0 

RAO-57. 2R57795 
CON-1. /RAO 
PI-3.1415926536 
THOPI-TViP*PI 
GN-1.4076576E16 
PT-6076. 11549 
RE-20925741. 
RnTATE-7.29211E-5 

ENTRY INITl 

ICNT-1 

SIGHA-l. 

XK-12345678. 

XHU-0. 

RETURN 

END 



O O uuu OCIUOUU 


subroutine input (TITISI, TITLES) 

THIS SUBROUTINE READS MOOIflEO NAMELIST FORMATTED DATA, 
IT READS A CARD ON UNIT 5» WRITES THE CARO Oh UNIT 6» 
WRITES THE CARD ON UNIT B (FIRST 72 CHARACTERS ONLY). 
THE TITLEl AND TITLE2 CAROS ARE NOT WRITTEN Oh UNIT fl 
BUT THE DATA IS PLACED IN THE TITLEl AND TITLES ARRAYS 
FOR TRANSFER BACK TO THE CALLING PROGRAM. TITLEl* AND 
TITLES* CARDS MUST BEGIN IN COLUMN 2 WITH NO SPACES. 

THE CALLING PROGRAM MUST BLANK THE TITLE ARRAYS# 

CALL INPUT(ARG1#AR62»)# AND REA0(8# INPUTO ) . NAMELIST 
DATA MUST BEGIN WITH tINPUTO AND END WITH tENO* BOTH 
BEGINNING IN COLUMN 2. 


C 

c 


c 


c 


c 


c 


c 

c 

c 


c 


c 


IMPLICIT INTEGER(A-Z) 

DIMENSION TITLEK8)# TITLE2(8)# CARD(8) 
DATA BLANK/iOH / 


REWIND B 
WRITE(6#A0) 

10 CONTINUE 

READ(5»50) CARO 
IF (EOF(5).NE.O) STOP 
WRITE(6*60) CARD 

BLANK COLUMNS 9 AND 10 
ENCODE (10#70»W0R0) CARO( 1 )# BLANK 
IF (WORO.NE.lOH TITLEl* ) GOTO 2U 
CARO READ IS A TITLE CARO 
ENCODE (72#8U#TITLE1) CARO 
GOTO 10 

20 CONTINUE 

IF (WORO.NE.IUH TITLES* ) GOTO 3U 
CARO READ IS A TITLE CARO 
ENCODE (72#80. TITLES) CARO 
GOTO 10 

30 CONTINUE 

CARO read is NOT A TITLE CARO 
BLANK COLUMNS 73-80 OF DATA CARO 
ENCODE (1C»90#CARD(8 ) ) CARO (8 )# BLANK 
WRITE(8*50) CARD 

IF (CARD(l).NE.lOH SEND ) GOTO 10 


REWIND 8 
CALL OATE(OAT) 

CALL TIME(TIM) 
WRITE(6#100) DAT»TIM 
RETURN 


^0 format 
50 FORMAT 
60 FORMAT 
70 FORMAT 
80 FORMAT 
90 FORMAT 
100 RQRMAT 
END 


(IMI) 

(8A1C) 

( 10X#8A10) 

( A8. AS) 

(R2»7A10) 

(A2*A8) 

( /////»10X#*DATE IS ♦A9/10X,*TIME IS *A9) 
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SUBROUTINE INTER (X,Y»NUM,8» A ) 

SECOND ORDER INTERPOLATOR 

SELECT FOUR DATA POINTS CLOSEST TO X TO INTERPOLATE FOR Y. 

X-INDEPENOENT VARIABLE VALUE 
Y-RESULTIN6 DEPENDENT VARIABLE VALUE 
LNT-Nn. OF ELENENTS IN A AND B 
B-ARRAY OF DEPENDENT VARIABLES 
A-ARRAY OF INDEPENDENT VARIABLES 
DIMENSION A(1S)»3(15) 

I-l 

IF (NUN.EO.A) GOTO 30 
IF (NUN, IT. A) WRITE (6»A0) NUH 
IF (NUN.LT.A) STOP 

C • 

IF (X.LT.AO)) I-l 
IF (X.GT .A(NUH-Z) ) I-NUH-3 
IF (X.LT.AO) .OR. X.GT.A(NUH-2) ) GOTO 30 
C 

LMT-NUM-2 
DO 10 K-A.LHT 

IF (X.LT.A(K)) GOTO 20 
10 CONTINUE 
20 CONTINUE 
I-K-2 

c- 

30 CONTINUE 
XC-A(I) 

Xl-4( lAl) 

X2-A(I^2) 

X3-A(I+3) 

Y11-((X1-X)AB(I)-{X0-X)*B(I+1))/(X1-X0) 

Y21-( (X2-X)*B(I)-(X0-X)AB( 1+2) )/(X2-X0) 

Y31-((X3-X)*B(I )-{XO-X ) *8 ( I'«-3) ) / ( X3-XC) 

Y22-{ (X2-X)PY11-(X1-X)PY21)/(X2-X1) 

Y32-C (X3-X)*Y11-(X1-X)*Y31)/ (X3-X1) 

Y-( (X3-X)*Y22-(X2-X)*Y32)/(X3-X2) 

RETURN 

C 

AO FORMAT ( //10X»*SUBRQUTINE INTER - VALUES IN INTERPOLATION TABLE •* 
• 13* MUST BE .GE. A*) , 

END 
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c 

c 


c 

c 


c 


c 


SUBROUTINE NEWTON 

ITERATES ON ECCENTRIC ANOMALY GIVEN A TIME (TCALU 

COMMON /CENT/ CETAf SETA»CG*SG>R*V»SEA 
COMMON /CONIC/ STAT£(6)» TCALL* KERR# STEP# T» T7ERO 
COMMON /MAN/ A • Al EE# AVG* E # OOG»E AR EE # ET A . I C IRCL * LOUS Y» P# 
A PGME #RREF#TANOM»TRUE* TREE #XN A, PER I on# EP 


KERR -.O 
K • -30 

OTH • TCALL - TREE 

EA > EAREE * XNA * OTM / RREE 

VALUE • AVG ♦ TCALL 

10 ECN • ( EA - E*SIN(EA) - VALUE )/(l. - E*COS(EA) ) 
EA ■ EA - ECN 
K ■ K •# 1 

IE (K.GE.O) GOTO 30 
IE (ABS(ECN)-l.E-9 .LE. 0.) GOTO 50 
GOTO 10 

30 WRITE (6#A0) 

AO FORMAT! / 6X 22H30 ITERATIONS EOR BODY // ) 

KERR-l 
GOTO 60 

50 EAREE - EA 

SEA > SIN( EA ) 

CEA > CnS( EA ) 

RREE • A ♦ ( 1. - E ♦ CEA ) 

TREE • TCALL 

STA ■ SEA * AIEE / RREE 

CTA > ( CEA - E ) * A / RREE 

ETA ■ ATAN2( STA. CTA » - TANOM 

TRUE • STA 

TRUE STORES STA EOR SIGN ON GAMMA 

60 CONTINUE 
RETURN 
' END 
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SUBPnUTiNE NOROAN 

GENERATES A ntBUALLY DISTRIBUTED RANDOM NUMBER. 

REAL NG.KX.MU 

common /BLRA/ MU# SIGMA# GAUSS# NGI16)# KX# ZRAN# ICNT 

TEMP ■ 0. 

DO IOC I-l#12 

KX • N6(ICNT)*KX 

KTEM3 • INT (KX/ltOOOCOCO.OO) 

KX • ABSIKX - FLOAT(KTEMPJ*100000001.00) 

ZRAN > KX*9.<jOOOOC 01CC 
ICNT • ICNT ♦ I 
IF (ICNT. 6T. 16) ICNT • 1 
iOO TEMP ■ TEMP ♦ ZRAN 

GAUSS ■ SIGMA*(TEMP - 6.0) ♦ MU 
RETURN 

end 
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SUBROUTIME TSTEP 

PROPAGATES ALONG A CONIC BY TINE INCREHENT <STE»*) 
SUBROUTINE CONIC MUST BE CALLED TO INITIALIZE 


COMMON /CENT/ CcTA,SETA,C6»S6»R>V»SEA 

COMMON /CONIC/ STATE(6)# TCALL* KERR* STEP* T* TZERO 

COMMON /CORO/ 0(3*3) 

common /dig/ ZEP0*0N£*TW0*RAD*PI*TW0PI*GM*FT*RE*C0N*R0TATE 
COMMON /MAN/ A*A1E6.AVG*E*OOG*EAREF*ETA*ICIRCL*LOUSY*P* 

A PGHe#RREF»TANOM*TRUE*TBEF*YNA»PERIOO,EP 

DIMENSION XU(3)* XV(3)* ZN(3)* 0N(3)* 0A(3) 

DATA ZN / 0«* 0.* 1« / 


IF (LOUSY. 6T.0) GOTO 80 
IF (ICIRCL.LE.O) GOTO 10 

CIRCULAR ORBIT REQUIRES SPECIAL CARE. 
eta • AVG ♦ STEP ♦ ETA 
TRUE • ZERO 
EAREF-ETA 
GOTO 20 

10 TCALL • TREF ♦ STEP 
CALL NEWTON 

IF (KERR.NE.O) GOTO 9C 
20 T« T ♦ STEP 


CALCULATE POSITION AND INERTIAL VELOCITY COMPONENTS. 

SETA • SIN( ETA ) 

CETA • COS( FTA ) 

R • HREF 

V - SORT( GMA(TWO/P - ONE/A) ) 

CG • DOG / ( R * V ) 

S6 - SORT( APS( ONE - CG * CG ) ) 

IF (TRUE.LT.O.) SG ■ -SG 

RCOMP ■ V ♦ ( CETA * SG - SETA * CG ) 

HCOMP « V ♦ ( SETA ♦ SG ♦ CETA ♦ CG ) 

CONVERT INTO INERTIAL X-Y-Z AXES. 

X - R « ( CETA * 0(1*1) ♦ SETA * 0(1*2) ) 

Y • R ♦ ( CETA ♦ 0(2*1) ♦ SETA ♦ 0(2*2) ) 

Z > R * ( CETA * 0(3*1) * SETA * 0(3*2) ) 

OXOT ■ RCOMP * 0(1*1) * HCOMP * 0(1*2) 

DYOT • RCOMP * 0(2*1) ♦ HCOMP ♦ 0(2*2) 

DZDT ■ RCOMP • 0(3*1) ♦ HCOMP • 0(3*2) 


XU(1) • X / R 
XU(2) - Y / R 
XU(3) - Z / R 
XV(1) ■ OXOT 
XV(2) ■ DYDT 
XV(3) ■ DZDT 
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nA( 1) «XU(2)*ZN{3)-XU( 3)*ZN(2» 
0A(24«XU(3»*ZN(1)-XU(1»*ZN{31 
0A(3)«XUa)»ZN(2)-XU(2)*ZN(l) 

AA-DA (1)*DA(1)^0A(2)«0A (2)-*'0A(3)*0A(3) 

IF (AA.GT.l.E-60) GOTO AO 
OA(l)-l. 

0A(2)-0. 

0A(3)-0. 

GOTO 50 
AO AA>S0RT(A'A) 

0A(1)-0A(1)/AA 
0A(2)>0A(2)/AA 
0A(3)-0A(3)/AA 
50 CONTINUE 

0N(l)-0A{2)AXU(3l-0A(3)*XU(2l 
DN<2» •0A(3)*XU<1)-0A(1)*XU«3) 
0N(3)»DA(1)AXU(2)-0A(2)*XU(1) 
AA«ON{ll*ONm+ON(2)AOH(2J*OH{3)*ON(3) 

IF ( AA.GT.1.E-6C) GOTO 60 
DN(1)-1. 

OK (2) >0. 

ON (3) >3. 

GOTO 70 
60 AA>SORT(AA) 

0N{1)«0N(1)/AA 
0N(2) -ON (2)/ AA 
0N(3»-0N(3)/AA 
70 CONTINUE 

VSORTH- XV( 1)*0N( 1) +XV( 2)*0N(2)+XV(3)*0N(3) 
VEAST»-(XV(l)*DA(l)+XV(2»*0A(2»*XV(3)*DA(3)) 

C U»OATE THE state TRAJECTORY VARIABLES 

H • R- RE 
GAM«PA0*ASIN(SG) 

7ET - RAO • ATAN2( VEAST# VNORTH » 

IF (ZET.LT.O.) ZET-ZET+360. 

ALAT • RAD ♦ ASIN( XU(3) J 

ALCN • RAO * ATAN2( XU(2)» XU(1) » - ROTA TE* ( T-TZ6R0 ) ♦RAO 
STATE (I )-H 
STATE (2) -V 
STATE(3»«GAN 
STATE(A)-ZET 
STATE (5) -ALAT 
STATE(6)-ALPN 
C 

eO RETURN 
C 

RO WRIT£(6»100) STFP 

IQO FORMAT( U lOX 16HTR0UBLE IN TSTEP AX 6HSTEP • 1PE15.5) 

return 

END 
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SUBPOuTINE uniran 
GEN tPATES" A UNIFORM «ANOOf< NUMBER 

REAL NG,KX»HU 

COMMON /BLKA/ MU, S I GMA , C- AL'SS » NG ( 16 ) # KX , ZR AN , ICNT 

IF (ICNT.LE.& .OR, ICNT.GT.16) ICNT-l 

RX • NG(ICNT)*«X 

ktEMP - INT (KX/iCOOuCOOC.OO) 

KX • A8S(KX - FL0AT(KTEMP)*10C000001.C0) 

ZRAN - KX*.nCCiOl3ulOO 

ICNT ■ ICNT ♦ I 

IF (ICNT. GT. 16) ICNT ■ 1 

RETURN 

END 
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RESOURCE REQUIREMENTS 
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Core 65k (octal) Tape none Plot no Graphics none 

DESCRIPTION 

Purpose: To calculate error statistics in the trajectory parameters, including 

an error covariance matrix, from flight experience. The covariance 
matrix is propagated on option. 

Input: Flight experience errors in the trajectory parameters. 

Output: Error statistics in the trajectory parameters, error covariance 

matrix, propagated covariance matrices, spin-stabilized stage boost 
covariance matrix, sensitivity matrix which represents the final 
covariance matrix eUEter propagation. 


Functional 

Description: 

Flight samples are combined statistically to obtain an error 
covariance matrix. The covariance matrix Is propagated using a 
Monte Carlo technique of sanq>ling, propagating the deviated 
trajectory along the conic, and statistically combining the re- 
sulting errors into another covarlemce matrix. The sensitivity 
matrix is developed from the eigenvalues and eigenvectors of the 
covariance matrix. 
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